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ELECTRICAL ANALOG COMPUTER FOR L LIMIT DESIGN N 


this | paper an analog computer capatile of ‘minimizing a 
function subject to large numbers of linear inequalities: is described. 
This” is the situation customarily encountered 
problem. 


tomary n- -space geometric c analogy. 


The problem of minimizing a linear function of many variables, subject al 


not impossible, to solve analytically for the number of 
associated with many typical realistic problems . For example, the large 
number of possible kinematic mechanisms makes the limit design of a typical — 
10- -member frame virtually impossible without the use of modern high-speed — 

id At the present it time | linear programing techniques exist for digital com- 
puters, but because of cost and the general unavailability, these machines can 


& used only for the most important design problems. 
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ntional computer lead toa machine 


It is the of this paper to a special-purpose machine, made 
up largely from available and relatively inexpensive components, capable of 
solving the limit design problem and, as a special case, the linear program- 

ming problem mentioned woe: a particular, this paper will be concerned 

with the two-dimensional, © redundant, frame structure composed of a finite 

umber of constant cross- section beams. The : area and shape of each cross- 

section may vary from beam to beam , subject to the restrictions implied . 

design problem is 3 stated i in the following manner: Given the general iis) 


the total weight of the structure, subject to the condition ‘that ‘the i ; 


“Limit design” “theory will be used in following and: an 
effort will be made to keep this paper self- contained, some of the implications — 
of wot basic assumptions a as well as s alternate approaches and solution tech- : 


assumed that structure is ot 
plastic members which deform entirely by bending action. 


(2) The linear density of every member in the structure is assumed to rl ts 


representable by the ne Saat function of of its s fully plastic moment nt Mj, i 
‘the range of interest. 


Based on ‘the assumptions, the total weight 


int 
- the see of minimization of weight is equivalent to the minimization of 
mechanism by the c combination of rigid bars and associated 
_ ideal “hinges” that can produce continued deformation without increase in ex-— 
; terior loads; it follows from the basic theorems of limit analysis that any 
-mechanism which everywhere satisfies ‘equilibrium, | as well as the flow law 
3 Dines, L. L. “System of Linear Inequalities,” Ann. Math. (2) 20, 191 (1918-1918). 
Jacques, Heyman ‘Plastic Design of Beams ‘and Plane Frames for Minimum 
terial Consumption,” Qu. Appl. Math. Vol. VIII, #4, pp. 373-381 asta, Pte 
Rg W. S. Dorn, H. J. Greenberg “Linear Programming and Plastic Limit Analysis - 
. e- Foulkes “Minimum Weight Design and the Theory of Plastic Collapse” Qu. Appl. 
Math. Vol. X, #4, pp. 347-358 (1953), 
ay: 7 J. Foulkes “The Minimum Weight Design of Structural Frames,” Royal Soc. Proc., 
ain ae R. K. Livesley “The Automatic Design of Structural Frame,” he 
i Mech. and Appl. Math., Vol. IX, Part 3, pp. 257-278 siateall 
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and its associated yield criterion, will lead to an admissible set of “unknown 
moments Mo} . The yield criterion is automatically satisfied by making the 
internal nada everywhere equal to or smaller than the fully plastic mo- 
ment. The flow law is satisfied by taking the as: assumed gee to have the — 
developed plastic moment of the proper sign. 
general, , most ‘structures will have an infinite number of admissible 
moment combinations. It is the problem of limit design to select that admis-— 
im sible set, or sets, which in addition minimizes the total structural weight. = 


the limit | load for a structure is of the elastic 


modulus of the structure, 10 


tural Contiguention for a given oak of loads is also independent of the elastic >. 
modulus. _ Thus, for convenience only, the rigid, ideally plastic structure will 
be used in the ensuing ‘discussion. ALL results apply, without change, to the h 
@® Instead of evaluating all unknown moments in terms of the revundants 3 
eras the use of the equilibrium equations: as was done by R. K. -Livesley,® 


\ terms of the unknown fully plastic moments as was done by J. Foulkes.6,7 in 
In writing these it will be assumed that the non- -dimen- 


‘The inequalities (see below) with such will be desig- 
nated as a group. In general, there are such groupings possible. The 
lowing statement applies to any generic group: 
ee for the system | of n _ unknown moments, any n of the virtual-work ex- 
pressions are taken as. ‘equations and the rest. as inequalities, with» internal 
s being greater than or equal to external work, the conditions of equi- iy: 
= librium and the flow law with itsassociated yieldcriterion will beautomatical- _ 
oT ly satisfied. "4 In the above, requiring that the internal work is greater | than the 
external work insures that the associated mechanism cannot form for the ad-— a 
missible set of unknown moments Mo; and presc ribed external tractions. 
Based on the three preceding sections, the analog described herein will 
minimize omg function 2 ‘subject to the virtual work inequalities 


Geometric Analog.— Expressed 
is viewed as a hyperplane in “n” is moment space. ” ‘Then the “ permissible” re- 
(no inequalities violated) be bounded by “k” ‘of these hyperplanes, 


bea 10 M. Zaid “On the ae Safety Factor for Plastic- -Rigid and Elastic- Plastic Bodi | 
es” un- 
a R. K, Livesisy “The Automatic Design of Structural Frames,” The Qu. Journal c of 
and Appl. Math., Vol. IX, Part 3, pp. 257-278 (1956). 


6 J, Foulkes “Minimum Weight the Theory 0 of Plastic 
‘Math. » Vol. X, #4, pp. 347-358 (1953), 
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it can be shown§ ‘that his permissible region is a convex set in the first t quad- 
rant. Representing function Q by its hyperplane | in moment 


self toward the origin. 
tangent to the permissible. region goa the permissible region is convex, | the - 
-tangency must occur at a corner. vl Included in 1 this tangency « condition is the | 
special case of the minimizing function being tangent to the convex region at 
‘Thus if a _ technique is evolved for examining the boundaries of the per- 
‘missible region, one need only compare the minimizing ‘function at each cor- 
ner of this region to obtain the solution, 
Since all groups have common virtual work equations at their common F 
boundaries, they will be joined together without moment discontinuities; i. e., 
_ if group A differs from B only in that m, > mg in A and mg > m, in B, then : 
ry = edge formed at mc = mg marks a boundary between A and B. "Thus me the 


in the following manner: __ 
‘— Q is fixed at a permissible value (Q), and, in addition, any n-1 of the 


of the n-1 equalities form an | edge, and the Q hyperplane they form a 


corner. A subgroup comprising the n-1 equalities and tthe violated inequality, 
“now > treated as an equality, is formed and, subject to the Qo F plane, all n are 
& automatically sampled and from then a new sub- -group of equalities are 
4 > that does n not ot violate | the remaining in inequality, and has a corner in the 5 ; 
entire group of k-(n-1) hyperplanes is again examined for violation 
_ with respect to the new corner (n-1 hyperplanes, and fixed Q, plane). Whena_ 
3 5 violation is found, the associated hyperplane is considered in conjunction with 
the (n- -1) hyperplanes and a new corner is chosen on the Qo plane that oes 
not violate the remaining hyperplanes, and lies in the original group. = 
33 Ir, _ By this process various corners are examined on the ~e until a per- 
at ‘missible one is found 


of the equipment of this | paper | r the above | process is ; automatic. widcnccaael : 

ie At the present time the convergence of this process has been proven for 

— only some of the general classes of cases. However, for all the ‘particular — 7 


nverge d. In the event that a ‘case should occur for 


3 Soe for a permissible corner, ‘the following procedure could be taken. 
) mn 1) of the planes « of constant m, inward until they form 2 acorner at this pare 
_ ii) Decrease one of the mj so that a line in n-space is traced in the direc- 


Shad 6 J. Foulkes “Minimum Weight Design an and the Theory y of Plastic mae oma 
Math., Vol. X, #4, pp. (2963). : 
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ii) Next, using (n-2) of the m; planes , subject to n to movement along the Qo . 
and the newly | attained boundary surface, one of the m > m j Planes is again de- 
creased, of the first quad- 
sz is traced in n moment space. Since the permissible space is convex in 
= first quadrant, this line will intercept a new permissible surface. = = 
(iv) By a continual process of intercepting new surfaces and discarding the 
my planes, a corner of the permissible space is finally reached. In following - ; 
‘this procedure care must be taken to properly consider the groups and the — 
associated | permissible boundary surfaces as the m; variables change in 
(3) The minimizing function Q is now decreased and a path 2 along a per- 
; - missible edge is traced until a permissible corner is reached. _ This is indi- 
_ cated by the intersection of the edge witha new hyperplane. 2 —t 
4) Of all the paths emanating from ‘the permissible corner, the computer 4 
automatically searches for one, represented by n-1 of the intersected hyper- 
4 planes, which with continued decrease of the minimizing function Q does not 
violate the nthhyperplane, 
(6) ) This process of moving along permissible edges from corner to comer 
any is continued until a corner is reached from which there are no permissible 
a ‘” edges in the direction of diminishing Q, or until a new group is reached. In * 
_ the former case the (Q plane is just tangent | to ‘the « corner and the solution i is 7 


_ (6) In the event that a group boundary is reached, | continued decrease of Q 
ts attempted along a path into the new group without ary mj | discontinuity be- ; 
tween the two ee. This proceeds in the same manner as in the former — 
7 (7) If further Q decrease is not possible into the new group, then the path a 
is constrained to move along the cc common boundary between the two groups. ba 
a (8) The above process is continued until there is no further direction of - ads: 
decreasing Q either along the common boundary of several groups or into a_ | 
_ new group. . In this case, ‘because of the convexity of the permissible region, 
the Q plane is tangent to a corner and a solution is obtained. r: a | ae ; 
_ It should be noted that a non-geometric, one-dimensional interpretation of 
a similar but non-automatic analog computer has been established by Ryder. 11 
In his discussion which was limited to the problem of a single group, Ryder — 
envisions two arrays of rods: positioned by the right-hand side of the 
equalities and the other by the left-hand side. Oneof the left-hand rods, which 
the Q function, is as far as which time the 


Basic: Components. this section the basic analog components will be 
ribed in terms of the operation that they are to perform. Their actual use 


i 11 Frederick L. Ryder “Linear Algebraic Computation iby Multi- -Winding Transform- 


Vol. 259, No. 5, pp. 427-439 (1955). 
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problem, but since it does not readily yield any additional informationandin 
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in the eomauter is left to the section on the “Complete Analog,” which which de- 
scribes how the actual computer, using only the basic components, solves the | 
limit design problem. It might be well for the reader to scan this section de- 
scribing components and then reread it in detail in conjunction with the de- . 
scription of the complete computer when it it can be’ better understood. In order — 
to keep the ideas of this paper as fundamental as possible the engineering de- 
sign problems have been avoided by indicating the components schematically. | 
It is left to the interested reader to fill in the engineering details of the actual _ 
equipment. In some cases this ' will result in the combining of several sepa- 
_ rate functions into one piece of equipment so that the final computer has even 
_ greater — simplicity than appears in this paper. In addition, in order to de- 
crease the equipment cost, it is possible to replace certain of the automatic 
by manual with some loss ‘of operating convenience | 


(5a) 


= Aj; 
‘The Ajj are prescribed in and the m,; are the non- un- 


known moments. The actual component used in this analog isa bank of multi- s 
winding essentially ideal transformers.11,12 This is shown schematically in 


q Fig. 1, and will be iiiiialis represented as shown in Fig. 2. The mj, 


1.-CUMULATIVE ADDER 
"winding. of of ‘these transformers, when Q and 
» bs linearly independent Ij are impressed and all other quantities are left 
free, there will be sufficient constraint so that all the m,; and remaining I, 
will take definite and consistent values, and negligible. current will flow. When * 
:: n-1- -s of the Ij are prescribed, the system is under-constrained, and ideally as 
Frederick L. Ryder Linear Algebraic Computation by Multi- ti- Winding Transform- 


_vol. 259, No. 5, pp. 427-439 (1955), 


2R.R. M. Mallock « ‘An Electrical Soc. 
A, Vol. 140, pp. 457-483 (1933). 
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"there is an s-fold infinity of eile values for the mj and remaining I,, with © 

negligible current flow. Actually only one set of values would exist but they 

- would be unpredictable | without a detailed knowledge of the very small trans- = 

former imperfections. When n-1+s of the Ij are. prescribed, the system 
over-constrained and in general equation (5) will not be satisfied. In addition a 
tend 
@ 
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(2) Element—(Fig. 3). “This uni unit of an amplifier 
that compares two voltages and, depending on the algebraic wo of the differ- - 


; that in - the ‘event that the amplifier is 5 disconnected, itr remains ; at its last | po- 
_ sition until it receives a new signal from te amplifier to change that posi- 


INPUTS 


Comparator a Energizer. 4). unit replaces 
large: number of amplifiers in a bank of constraining elements. This replace- = 
is accomplished by the use of stepping switches (shown schematically | 


4 pair of rotating contacts) together with a reversing switch to feed a single _ 

Since there isa a requirement to open as well as close the relays, and since 
amplifier designed to apply only a large pulse for a small voltage 


ali 
— 
— 
| 
a. 4 


a reversing switch in order to discriminate between an opening and closing 
oo. In addition, since the relays are of a two- -~separate-coil type (with — 
’ memory), one coil being energized by the apeuing signal and the other by the 


IF"O" AND Ij < 


"n n+! INDICATOR" 


closing signal, it is necessary to introduce an open- 1-close relay to ‘iain the 
impulse from the amplifier to the correct relay coil. This open-close relay | 


the reversing obtain their signals fr from the “ ‘nor n+l” indicator, to 


FROM INEQUALITY | 


+ Maintainer. —(Fig. This unit supplies the relays ina bank 


_ of constraining elements. The amplifier system supplying the impulse to the 
"relay coils is that described above in section (3) (Inequality Comparten. 


— ined polarity, it has become 
— 
— 
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LIMIT DESIGN 


MIT DESIGN 
Because of the fact that each relay receives only a single pulse to open or 


close and there is no signal until that relay | is reconsidered, | 
latching relays with two coils. 
a4 


it is necessary y to ; 

The function of the relay is to impress voltage Cj if the associated Ij ever a 
tends to ) violate inequality (4); otherwise, no such constraint is applied. The 

in the I; line maintains slightly different from C;, j, so that even 
_ when the line is closed a decision on whether or not to keep Cy impressed i &§ 

Ss be made on the basis of the polarity of the voltage drop across this | resistor. q 


In addition, it is bie have whether or not any particu- 
lar isimpressed. 


uw 
it should be noted here that in a design it is to compare > 


: -C; with zero (ground) and by s so 9 doing eliminate one set of contacts in the 
inequality comparator (3), 


(5) “n or n+l” Indicator.—(Fig. 6). This unit is used to indicate whether 


the system is exactly constrained (n equations impressed) or 
fa + 1) equation impressed). 


that the potential between one B is zero when n-1/2 resistors are (theo- - 
-Tetically) by relays of the maintainer. (Since the 


— 


EQUALITY MAINTAINER 


FIG. 6. OR n+1 INDICATOR 
equation is always maintained this corresponds to ) n+1/2 equations). The 1 am- ii 
as plifier is designed so that for a potential difference A-B <0 there is no out- 7 
ig 


put, and for A-b =: 0 there is an output. _ Thus for n-1 or less relays (n equa- 
tions) engaged there is no output and for n 


(n+1 equations) or greater there is es: 
an output. This signal is fed into the relay of the inequality comparator to es- 


tablish the opening or closing and reversing circuits as described. ; fe 
oe (6) mj Limiter. —(Fig. 7). This is a relay bank that used. the inequality — 
- comparator in order to set a lower limit on m; . The lower limit is shown as ¥ 
_avariable source, which is introduced to the inequality comparator together 
with mj in the same manner as C; and I, . In addition, by means of patch cords 
any of the m’s can be set equal to each other if required. As in the equality =e 
maintainer, a light indicates when the lower bound on mj is in operation. ere. : 
Ty (7) Group Indicator and Energizer. —(Fig. 8). This is abank of constraining — 
elements: connected in series so that any ordering of the quantities can be 
‘specified. the actual ordering is as specified, then a voltage is impressed 
4 on the ith output line . Fig. 8b shows a typical arrangement of the wiring or 
group i; ; the ith 


switch thrown when the i group 
switch it in Fig. , 8a | ‘is set to the “ on” or* over- ride” position. iu 


4 

_ 

— ae 
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This unit has an over 


- when the n or n+l indicator gives a signal indicating over- ~determination he 
over-ride automatically functions and the previous groups rer remain in energized. — 


INE equaity | 
COMPARATOR AND 
ENERGIZER 


FROM INEQUALITY 
COMPARATOR AND 


be 


5 ee the u unit shown, more than one group (ordering) « can 


any one time. which group or groups are function- 


mee 


_GR 


a —(Fig. unit is a of multi- 
receive their from the group indicator and energizer, any par- 
ticular group activating only its assigned relay. The relay acts asa switch, 
its so as to take the inputs and the Cig 


— 
— 
— 
— ate — 
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ze selector can be thought of as supplying outputs C; or zero voltage as_ 
required by the group being considered. _ It should be noted that if the ampli- © 
fier design required it, , Zero voltage could be assured by grounding the proper 
pole on the iia instead of allowing it to ont. . 


m > > 


FIG. 8b.—TYPICAL RELAY AND AMPLIFIER WIRING FOR 
GROUP INDICATOR AND ENERGIZER 


(10) Re: Reader. .der.—This is a unit that has one or ‘Several voltage- reading 


GROUPE GROUP GROUP 3 


FIG. 9. —GROUP SELECTOR 
‘:Oneatinie renee —The schematic for the complete analog using the vari- 
ous units described is shown in Fig. 10. The operation will now be explained 
_ in terms ofthe theory discussed above. = 
‘The Limit- (Eq. 4) are first wri itten and 


— 
— 7 _ each of which is either the voltage C; or zero current (obtained by opening the | i ae 
line). ~The amplifiers using this output are so designed that they respond to 
— 
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laxing this to establish a permissible corner. 

a permissible | set of “n-1” m | values | is set on the mj limiter, using the _ 

limiter over-ride switch to these settings. An associated permissi- 
ble Qo is set on the Q generator, so that the groupand n of the m; are defined. — 

7 = The proper group switch is set on “over-ride” in the group indicator and 
energizer, and the group selector automatically sets Cig equal to ) either Ci or 
zero depending on whether or not the inequality is in the group to be —— 

sidered. _ The inequality comparator finds no violation since the mj are per-— 
missible. At this time the m; limiter over-ride switches are disconnected; 
which brings 1 no further change in the analog. 
, _ Now, assuming that a group is maintained, one of the m; is decreased (its { 
oe plane is moved) until one of the inequalities just | becomes violated i (the line in 
space just touches a permissible surface). The inequality ‘comparator 
energizes the proper relay in the equality maintainer and the associated Cjis — 
engaged. This action is automatically relayed to the n or n+1 indicator, which > 
starts the opening cycle in the inequality comparator. Now it has been indi- 
 cated11 that at least one of the n-1 over-determined quantities will supply the 
inequality comparator amplifier with a potential difference of proper = | 
to open its relay. Thus this value of mj is automatically discarded. _ This im- > a 
Ss that any further movement along any of the lines emanating from the | 
corner not containing the discarded plane will not violate this plane. The dis- aa 
carded mj generator can now be set at its lower limit value. In this same _ 
manner the m; generators can be decreased and eliminated one at a time until 
all of the original “n-1" 1” constraints on the m,; are replaced by equalities — 
At this time a corner thas been reached on the permissible surface. 
= he Q generator is now decreased (the line subject to n-1 planes deter- — 
"mined by Ij-Cj=0 is followed) until the permissible surface is intersected. 
- ‘This forms a new corner. The ‘system now becomes over-determined and, in 
= same manner as described above for the discarding of the m; plane, one | 
of the non-permissible I;-C;=0 planes is automatically eliminated. Further 
‘decrease of the Q generator (plane) will not violate the discarded plane. This 
- process of decreasing the Q plane is carried | on, and the analog automatically. , 
selects and discards ag] so as to travel from corner to corner on the per- ie 
-missible surface. Finally, a value of Q is just reached for which the system 
remains over-determined (the n+1 light remains on) or also just leaves the 7 


group. _In the first event, the minimum value of Q is attained and the problem | i 5 


iin the latter event (the group changes), in order to preserve | conditions that 


setting in the group indicator and energizer. ‘Thus, even though mace 4 is an 

indication that the group has changed, the group selector maintains the previ- 

_ ous group energized. At this time the m; ; lower limits are manually raised to 3 

the permissible values at the last corner of the group boundary to prevent any D, og 
_ accidental loss of this position. This is best accomplished by raising the low- - 


er limit values in the mj generator until the associated -Ci= Planes are 
11 Frederick L. Ryder ‘Linear Algebraic Computation by Multi - Transform- 
ers” Vol. 259, No. 5, pp. (427- 439 


set on the Cj generators. With these coefficient settings the problem is 
or brevity this section will contain only the discussion of the technique of 
ee 
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just maintained, as can be by observing the equality lights of on the « equ ual - 
Now all the groups are examined, one ata time their switches 
to the “on” position), until a new group indication is obtained. This process 
% ‘cannot over-determine the system because the non-applicable groups cannot — 
energize the group selector. In addition, the new group is entered in a con- 
tinuous manner so that the same common permissible equations are veered 
= in the new group at the same permissible corner. The switch for the pewees . 
group is now set to “off”, and the new group switch is set to“over-ride”, the | 
& lower limits are set, and the Q generator is diminished just as in the 
= the event that it becomes impossible to enter a new group group, (the n+1 light > 
os remains on), it is necessary to consider the groups formed by the combina- 
of the groups. . This is is accomplished by “making the 


‘FIG. 11. 


This pr process of entering new groups is continued until the aiid of ~ es 
reached at which further decrease leads to an over- This 


Given: Structure shown in in Fig. 11. 


Find: Values of m and m that iz 


‘It is first to set up the that will will be used in in the 


From Eq. 2, the weight function 


= my + 


4 
a4 
a 
4 
4. 
“ts 
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aot 


2 


1 + 2mg + 2mg > > 


e 


_ sO that (using superscripts to indicate group) 


— 


) A set of | pars values of m, = Mp = = 0.5 and Qo = = 1.75 was selected. This — a 


resulted in mg = = - 0.25 with no inequalities violated. The associated a 
group was the combination of groups 3 and 4, but no attempt was made . 


=1, 75 and ‘M2 = the value of my, was de- 
face in A. pesgesernr group. This resulted in m, being decreased to i 
0.25, and Mg correspondingly ‘becoming 0.50. The group path was as 
; follows: group 4 was exclusively entered immediately on reducing my 
and remained until a new group formed from groups 2 and 4 was en- bi i 
‘tered just as the ers surface was intersected. The intersection ? 
=1 in in group 2, and and 2m, + 0 3 aie 


(2 10) 
— ‘Using the a anually (basica = 
— only, 
= 


both 2 and 4 my = Me and seeing that 
mg > mj, and enforcing one of the intersected permissible planes (the 
i other melee identically satisfied), the Q generator was decreased. _ This 
process was continued until a new permissible plane in the group was 
“1.87, ‘my = Mg = mg = 
First 2 ‘and following the 
two associated permissible surfaces for the direction of further move- 
4 = it wae: found that the slightest decrease in Q led to an intersection 


ix 
Psi of the two. associated permissible surfaces, it was found that any de- 
- crease in Q led to an entrance into group 2. Thus it was a to + 

proceed into either of these groups exclusively, 


© Selecting the common | group formed by the combined groups 2 and a 


— oo its associated two permissible surfaces (there are actually two | 3 


4 


ie identical pairs of permissible surfaces), it was found impossible to de- 0 


a] om rease Q subject to o one of the e permissible surfaces and to the require- 7 
- ment of this new group that mo = mg, without violating the remaining 
_ permissible surface. Thus no further decrease of Q was possible, so 
_ that the values indicated above correspond | to the minimum- -weight § so- + 
lution . Itis of interest to note that, even with the manually operated 
machine used, the actual running time to go through the entire ae, 
_ as indicated in steps (1) to (5) above was of the order of 10 minutes. a 


3 ‘For the automatic machine described in this paper, the ‘sctual running 
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LIMIT DESIGN 
NOMENCLATURE 


= equals C; or 0 


= number of groups; 
= length of member; 


- fully plastic moment of jth member; : 


— 
m; non-dimensional fully plastic mo vent; 
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Three commonly used systems of shell equations are employed _. 


~ solve the case of hydrostatic pressure acting ona segmental arch 1 ring. re 

The x results with solution by classical approach anda 
classical approach ‘including certain modifications. . Serious limitations 


of Donnell’s equations are discovered and the nature of differences be- 
With the increasing interest in application of shell theories to the design 
of arch: dams as exemplified by a whole series of papers(1,2, 3,4,5 6)2 al- 
‘ready published and further work in progress in various parts of the world, 
it becomes important to examine which of the various systems of equations | 
a proposed is most accurate and what is the > order of differences that may” p 
arise. “Studies of this nature are by no means new. Within the last decade _ 
‘several articles discussing this matter have appeared.(7,8,9,10,11) In all of 
these however, the accent was on the comparison of the roots of the charac- _ 
teristic equation which ; is used in the solution of cylindrical shell roof prob- ed 4 
— Very little information is presented on actual numerical comparison _ 


of stresses and deflections developed. It was with this in mind that the present _ 


Note. —Discussion open until June 1, 1960. To extend the closing date one month, a 
ities request must be filed with the Executive Sec retary, ASCE. This paperis part | 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the : 
_ American Society of Civil Engineers, Vol. 86, No. EM 1, January, 1960. [ear ; 

~ Assoc. Prof. of Civ. Eng., Northwestern Univ., Evanston, Il. 


Numerals i in parenthesis —thus » (1)—refer to corresponding items in Bibliog zraphy. 
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cylindrical shell) subject to uniform radial pressure was solved and as a 
sult, in addition to the numerical comparison sought, some general 
conclusions have been 

doy As the problem, though simple, is precisely of a kind applicable | to the 7 

arch dam, a ratio of radius: to thickness considerably greater than occurring 


_ by classical methods of structural analysis. (12,13,14) A comparison a 
fore of the shell theories with the and some 


te 
a SHELL EQUATIONS AND THEIR SOLUTION 7 
‘Three commonly used types of equations are those attributed to Flugge, 
and Donnell respectively.(15,16 917) first are usually ad- 
a to be the most accurate, but are » seldom used in ‘their. complete form — 
3 designers. Timoshenko’s equations form the basis of the ASCE manual(20) © 
7 while Donnell’ S equations are ‘very popular with shell designers in Britain. _ ie 
‘Some comment, perhaps, on the assumptions in the derivation of 
‘ln equations may be made here without acting greatly to the length of _ 
present paper. 


ion the behaviour of of thin ‘shell elements: ss ar 


Plane sections, the he middle “surface, remain n plane during» the 
BF esxe= of deformation and normal to the middle surface. 


ii. The deformations are small and therefore their effect on the configura- 


tion of the element for which the equilibrium equations are 


‘to th the | middle plane negligible effect on 


ited of t the shell is small with its of 
_ While most of these assumptions are analogous to those usually employed | 
is conventional beam theory and i therefore subject, presumably, to the same 
: limitations, differences arise in the manner in which the last one is inter- 
_ preted. Without going into the full details of derivation it can be at once ~ 


recognized om, , in integration of the for the various | stress + 


resultants, terms ‘the type —— occur. In these 2 is the perpendicular 


esr 


— 
ua 
4 
— 
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limit the i 
‘referred to as D. J 


ARCH RING 
ze - of tl the correct value.4 4 On the hand Fligge expands a 
expressions neglecting terms of an order higher z2, Although the ap- 
_ proximation ; appears to be better as the neglected terms are small, it is pos- — 
7 = sible for the relative error to be large if the retained terms are themselves - 
small. Thus, on occasion, it may be possible for the bending moments given 
by the different theories to have widely different values if the moments are 
7 themselves of small magnitude. This will be apparent from the considerably — 
7 different expressions for these which will be quoted later. rr a, 
ae’ _ The approximations used in the derivation of Donnell’ s equations are not 
as consistent as those described above. In the various expansions only the © 


lowest significant terms | are retained and in addition some terms from the 
equilibrium equations are arbitrarily dropped. Although arguments are 
4 quoted for the logic of ‘such approximations, (9,10) their justification must lie | 
in the accuracy of the final result. it will be shown later that under some 
circumstances completely fallacious ‘results can be obtained by their use. 

- It is not proposed here to quote the full differential equations which govern 
the normal, tangential and a: axial _ displacements of a shell as they are some- 7 
what lengthy. ‘These equations can be easily found in the references quoted. 

_ In the special case of a shell barrel of great length, or alternatively ofan 
arch, ring, the equations reduce to a relatively simple form. If Poisson’ Ss 


‘pespectively then the ‘equations take the following form when 
normal loading p is acting: 


Timoshenko; (wooo O00) -w)=p 


Ne DI 


‘ A small err error o mated or iginally in Timoshenko’s book. ‘This 
is perpetuated in| the AS 20e equations 81, 84b and 84c/ In this the 


"equation 81 read the corresponding terms in the other 


equations should be amended. The error was out by Am. 


Tr 
— 
i ay ie 76 p 160 1954/ In the present paper the results are not affected by this error. = * 


| 


‘Fig. 1 defines the origin of coordinates as ditt as s the positive directions of w ae) 
and v. M and N represent the moment and axial thrust as defined in Fig. 1. | 
i a is the normal pressure reduced to the action on the middle surface of the 
ae - That is, if in fact the pressure is acting on the external: surface its 


value should be increased in the 2. any ofthe sub 


‘With: p la ‘simple. case of ‘all the 
above equations can be solved exactly by orthodox means. It can eedinn be — 
checked by the reader the following are’ the general solutions: 
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‘ARCH RING» 


pee. of t the pow of the arch supplies s six err conditions from which a 
- six constants can be evaluated in each | case. . The algebra involved is tedious © 


but not difficult and anes have been obtained for various numerical para- -— 


‘Timoshenko’ s ‘equations are shown in Table 


“FIG, 1. 


The solution 0 of Donnell’ s reduces to a fairly simple form. Rio 


| 
~ 2 ag) * 
12 (5g) + 12 ( 
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‘tome ee: the Equations .—An inspection of Donnell’s solution for 
the case of § = 7, that » is, the problem of a complete ring, shows immediately 
o that it cannot be correct. This completely symmetrical case should result in 
-- moment which does not vary ry with the position. The reason for this not being 
becomes clear on further examinations. Donnell’ equation is not satis- 
‘fied if a simple rigid body is "imposed on the whole shell. With 
a displacement shown in ‘Fig. 2: 


TABLE 


_ Substitution of this into equations (3) ) show; that when p = 0, that is a case of . 
no load, neither equation (3a) is satisfied nor is the moment equal to zero at | er 
all points. On the other hand both Timoshenko’s and Fligge’s equations are _ 


. are caused by this rigid body movement, a condition essential to a proper 
formulation of any elasticity 


all of which A= 1+ — ( 
A 
— 
— 
q 
identically satisfied by the above relationship and no internal stre 
i 
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= a Timoshenko -0.01721 .00613 0.00210 0.00717 ‘00888 
Mote Fligge -0. ) 0.00133 0.00640 00811 q 
PR Classical -0.01618 0.00232 00718 0.00882 
Classical 3 0.00149 0.00635 =—-0.00799 


Donnell -0 91000 .9100 -0.9100 -0. 9100 9100 
Timoshenko + -0. 9370 0.9259  -0.9174 9126 -0. 9109 
Fligge —-0.9370 -0.9259 -0 9174 -0 9126 -0.9109 
Classical  -0.9396 -0.9290 0.9211 0.9163. -0.9147 

Classical -0.9396  -0.9290 -0. 9211 -0 9163 -0.9147 


‘Donnell 96001 5000 1.7070 
=WpR2 Timoshenko 0.3244 0.9658 1.5209 1. 7347 


Timoshenko -0.1602 -0. 2158 -0. 


Timoshenko 0.00962 -0.00353 0.00115 0.00409 0.00509 
Mo = 2 Classical -0.00939 -0.00344. «0.00112 0.00399 0.00497 


-0.01022 -0 00427) 0.00029 0.00316 0. 00414 


Donnell 9698 -0.9698  -0.9698 -0.9698  -0.9698 
N' 


Timoshenko 9853 0 9792 - 0.9745 9716 -0.9706 


wee 0.00034 0.00327 0.00428 


Classical -0.9857 -0.9779 | 9751 -0.9723  -0.9713 


Fligge -0.9853 9792 -0.9745 9716 -0.9706 7 
= 


Donnell 0.3479 (1.0229 1.5981 1.8184 
Fligge 3393 1.0201 1.6158 1.8468 


Donnell -0.1983 --0.2975  -0.2027 


Timoshenko 0 (0. 2256 -0.3055 -0.2093 
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wets 


-1.0000 


Classical 


3 4 
00173 -0.00053 0. 0. 00852 
"900137 «0.00039 00160 
0.00225 - 0.00044 0.00077 
0.00038 0.00158 0.00200 


-0. 00223 ~0.00045 00075 00075 0.00117 


-0.9938 -0.9938_ -0.9938 -0.9938_ 
0.9961 0.9948 -0.9944 

1.0000 -0.9979_ -0.9961  -0.9948 0.9944 

1. “0.9979  -0.9961 0.9949 a 9945 
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RING) 

In the case in point just such a translation is imposed on on the whole — 


occur in the central positions, these equations will lead to serious - errors, 
‘this indeed is very clearly shown inapaper by T. Moe.(7) 


A Numerical Comparison. — —In Table I moments, ‘thrusts 2 and displacements. 


are given for the. case a= —= = 10 and different central angles. 7 It will be noted 


that while the normal thrusts in the arch ring do not vary appreciably this is 


not the case with the bending eae? -In ‘Figs. 3a to 3d these results are 
_ Comparing deflections, it will be ‘noted that the results as given by — 
‘Timoshenko’ s and Flugge’s theory are practically identical; Donnell’s ay 
however, gives some discrepancy which becomes larger as the central, angle 
of the arch increases. In none of the cases is the discrepancy very large, 


- reaching a ‘maximum of about 8% of the largest deflection in the case of . 


On the other hand the bending ‘moments , derived by the various theories 
differ to a very appreciable extent, especially for central angles greater than 


_ 180°. A glance at the figures: will show the large errors in local moment and 


‘It is of interest to note ‘that in the case of Timoshenko’ 

theory yields a sero bending moment while Flugge’s theory yields a —s 

_ negative moment. The explanations for this lies in the thick ring theory. 
is well known that in a case of a thick ring subject 7 
‘the distribution of stresses follows the well known formulae(21)_ 
which result in higher stresses on the intrados and lower on the extrados. 
Although in the actual evaluation of stresses the | usual linear bending moment — 
formula is assumed, the stresses calculated by Flugge’ s result agree very 


_ closely with those obtained by the thick cylinder formula. — In the above case > 


10) these ‘stresses are + ‘different from average s stress” (which 


would be obtained by Ss In this instance Flugge’ 

theory gives results nearest to the correct solution. The differences 

7 7 however only of the aii associated with think ‘ring theory, and will be e negligi- 
ble in thin shells. 


q tt may be argued that the relatively large @ifferences | in the bending» mo- 
ments given various theories result in only of stres- 


(hoop) stresses. In ‘the example computed here, the discrepancy between | 


as can be seen if the expressions are. examined; 
_ differences of stresses as oe md the various theories diminish, also. = 


cin 
> before and throws serious doubt as to the range of their applicability. 
“In the cage of long culindrical shella where large cumulative flectiona can 
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am stresses caused by | al uniform temperature change. As will be seen later, 

the stresses in this case e are essentially to moments calcu- 
general technique of stresses in will be dis- 


cussed elsewhere. (2 ) Application this to the problem in hand 


bp) As in fact both the inner and outer boundaries are unstressed » an ordi- 


oR? use of previous results 4] 

‘The results of edilibenar a of stresses ofa and bis shown in Figs 4a to a 


; AS before, it is of iatabest to examine the limiting case of a complete ring. 
a In this case Timoshenko’s solution results in no stresses but both Fligge’s 

and Donnell’s equations result in non- zero states of stress. This is an — 


strained body is not ‘capable o of producing any iy stresses and it appears ‘that the 

very same terms which previously made Flugge’s equation result in more> 

correct stresses now cause this error. This point | has recently been com- 

mented upon by Eisentrout(23) who corrects the basic equations to. take 
B account of the changes of thickness in the direction normal to the plane of the E 
- shell. In the case under discussion such changes clearly may be of impor- 
tance due to the internal and external opposing pressures appliedin b. 
Although it was "mentioned before that Fligge’s solution results ap- 
g proximating to ‘the thick cylinder theory, in this case the errors are large 

enough to cause a non hydrostatic stress distribution whea uniform pressures > 
4 are applied to the inside and outside faces of a thick Bis ye to the re- _ 

quirements of the Lamé theory. 
7 CLASSICAL THEORY 0 OF . ARCH RING 

which has been solved is obviously amenable to standard 

classical procedures (such as Castigliano os Vines work). In fact, because © 
_ of its importance in connection with the study of arch dams, numerous solu- 

; tions have been obtained and tabulated in the literature, as previously men-_ 

, tioned. Clearly in this approach to the problem the ‘contribution of axial 

| ao forces to the deformations is of paramount importance as these comprise the _ 

primary action of the arch. The relative importance of the shear distortion > 
is difficult to assess “a priori” but in the analysis which was performed this 
_ term was included. _ The effect of the shear was found however to be almost 
negligible in all results. The calculations which shall not be reproduced here _ 


a to the following expressions for on ‘moment and thrust (a - being the 
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Where K = 5/6 is is the shape factor used in in cal 


resulting moments and thrusts have been computed we. are 
cluded in Table ‘Figs. cy 


and 4 show that the agreement of the classical 
theory with Timoshenko’ s shell solution is almost exact; (the very small dif- 


ferences being due to the shear term effect). 
a: In the conventional _ “classical” theory of structures, it is assumed that the | 
rotation between adjacent ‘sections is caused purely by the action of the mo-_ 
“ment; and the axial elongations of elements are due only to the action of the © 
Such a comparison was recently discussed L. (ref. 6) in which he 


EM1 
q 4 
4 culation of the shear distortion. 
a. i[ [ta 
— 


axial forms. That this 7 ‘only an can be seen 
-*Fig. 5. Even if plane sections are assumed it will be observed that the axial 


stress causes not owe an the plane c of magnitude 


(positive with the direction the moment). 


oe ly by the reciprocal theorem the moment M must in turn cause an elongation 


of the element of a an amount AER® 


That this the caseis not immediately evident, but can if Winkler’ 
_ theory of bending of thick sections is considered.(24) These additional de- 7 
_ formations can be coummares in the classical analysis(24) and an analysis of 


the example been conducted on this 


support Cro 


= 


expressions Mg and Ng are essentially similar to expressions (9) with the 
main difference lying ina constant term - = { ; 
+The corrected values of My are shown in Table I and are again plotted on 
Figs. 3 and 4 - It is seen that the answer here is almost identical with the re- 
sults obtained by Flugge’ s theory. of the essential differences between 
ss Timoshenko’s and Fligge’s equation now becomes apparent. It will be seen . 
that, if the derivation of both these equations is examined, Timoshenko ~ ol 


essentially the same simplifications as are conventionally made in sacral 
theory while the higher order terms retained by Flugge are to some aly 
- to considering the thick ring (Winkler’s) bending theory. It is F 


be expected t that it ‘in thin: shells the differences due to this become aan 


—_— © — 
— 
— 
ati the element of an amount - —— ds, 
also a relative rotation of the 

— 
7 
— 

4 

ete 
— 
Summari Mpoints should be reiterated: q 


theory. As such differences decrease with the | decrease of the 
t/R ratio, they are not likely to be of importance in reasonably thin 
‘Shells. On the other hand, in thick shells the neglect of shear distortion 


When computing thermal stresses, the use of Fligge’ relationship 


_ tends to lead to less accurate results then those obtained by the simpler — 
qd) Donnell’ suffers from the that rigid body 
_ displacements may cause appreciable stresses and should be used with a 
caution. In the case of arch dams, as the central angle is small, this — 
does not appear to. cause serious differences and in fact its relative 
‘tne may be a considerable advantage. Although results of cen- 
tral deflections of arches due to uniform water pressure are tabulated 


5.—ROTATIONS CAUSED BY AXIAL FORCES. 


in various publications mentioned, the ex extremely simple form which can 


derived from equation (7a) by putting 


crow 

prove useful in preliminary centilever” type of arch dam 

analysis. (it must be remembered here that p stands for the pressure 
transferred to the middle plane of the arch). This expression will be 


found to agree within a ofa a per cent with results of published 


e) Classical theories show very good agreement with shell a naming 
_ The author would like to thank Mr. A. L. Parme of P. ota sah sel 


his interest in the 2 problems discussed in this yer. 


erie a) Displacements calculated by the various shell theories show very good © cen 
b) The essential differences between Flugge’s and Timoshenko’s equations Pins, 
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RELAXATION THEORY OF CREEP OF METALS 
By Francis H. Ree, 1 " Taikyue Ree, 2 and Henry Eyring? Sina 


SYNOPSIS” 


ola A model for explaining ‘secondary creep at bia temperatures is |  pro- 
posed. Applying rate-process theory to this model, an equation for rsec- 
ondary creep rates is derived. This equation. is applied to creep of 
aluminum, aluminum base alloys, nickel with good results. 


is well creep consists of a primary secondary 


Note. .—Discussion open ‘until June 1960. To extend the closing date one month, 
Bere request must be filed with the Executive Secretary, ASCE. This paperis part _ 

f the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the 

American Society of Civil Engineers, Vol. 86, No. EM 1, January,1960. 
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_ Attempts to correlate the secondary creep with dislocation have 
een’ made by J. J. Kanter,4 Frederick Seitz and T. A. Read,5 Walter J. 
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January, 1960 j 
Kauzmann,§ , A. S. Nowick and E. S. Machlin,7 John E. Dorn, 8N. Mott 
ie Weertman. 1011, 12 12, 13, 14, 15 The five last-named investigators neve-ah- 


to apply the t ‘theory of rate processes16 to 
PROPOSED ‘THEORY 
= Model for Creep or Creep.—In order to explain creep behavior at high eneniivee, 
the activation heats for creep | for many metals equal those for self- 
diffusion,’ the following model is assumed. 
a ) A dislocation is generated at a Frand-Reed source (F-R) in Fig. 1: asa 
a result of the applied stress, f, and starts to glide on the slip plane, F- -R os, , 
which includes the dislocation and the slip vector. 
_(b) The dislocation is held up whenever certain bad sites lying in the éite- 
cation encounter other bad sites on the opposite side of the slip ge ll 
be be called “the counter plan” here. In order overcome this obstacle, 


oly reaches a sub- -grain (or grain) boundary, B with a resultant displacement, 

.. oa shown in Fig. 1. In Fig. 1 a plan view of the slip plane, F-R OB, is % 

_ show n, where t the dislocation line is represented by a line of d' e.. and the bad 

¥ “sites by crosses. The dislocation line, d'd', is located on the slip plane (the 
plane of the paper), the counter plane lying lower by an atomic distance. The 
bad site, which | belongs to the dislocation and is held ty? the bad site, a, s 

on by at. . Actually | the area shown in the plan, A , is the area over 


which the dislocation has slipped to reach grain with a 


In Fig. 2 (a) view of aa’ is shown. the nature of a bad site 
is not entirely clear at present, it is probable that a change in the charge “gst 
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“Theory of Steady State Creep Based on Dislocation Climb,” by J. _Weertman, 
11 “Creep of Aluminum Single Crystals,” by J. Weertman, J. App. Phys. 7. 
“Creep of Nickel »” by J. Weertman and P. Shahinan, “J. Metals, 


27, 1956, p. 1189. 
«14 “Steady-State Creep Through ‘Dislocation Climb,” by J 
Creep of Crystals,” by J. Weertman, App. Phys., fol. 28, 
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-eath other of the » bad sites. In any case, the held- -up concentrates 
_ ae = on the neighbors of the bad site, a’. . This stress is evidently relieved — 
a by diffusion of neighbors of the bad sites, since the ; activation heat for creep ie 


ig (d) Self- -diffusion involves two processes: the formation of a vacancy and > 


a hg the jump of an atom into the vacancy. In Fig. 2 (b)the free energy barrier for Bs? "4 
self-diffusion is indicated diagrammatically. (Actually free energy is re- 
quired and is defined for points ona surface. Here, 


‘a 


reaction is hindered as x Here, the pon free e energy for 
+ the forward reaction is reduced by w= A2 AZ {/2, and that for the backward 
c _ reaction is increased by the same amount, where the product, A2A3, is the -: 
fective area occupied by the diffusing atom. 
Rate of Creep.- —Since the gliding of a dislocation on a plane without 
cles is the creep rate, 8, is determined by the rate a 


— ie 
— 20 
— 
"4 art is the activation free energy for the pro nee = 
ug 


dislocation over bad sites idles in turn, i is determined by the rate of self- 
- diffusion. Let Z be the number of encounters between bad sites which the dis- a8 
_ location experiences while it travels from a Frank-Read source to a grain i 
Ea boundary (Fig. 1), and let k"' be the rate of relaxation of one of these en- ae 2 
counters. 


During the time of Z/k", one ‘dislocation through the he slip 


Reac tion “Coordinate 
F G.2 2.—(A) DISLOCATION HELD UP BY AN ENCOUNTER OF BAD SITES. 7 
the F-R source to the grain yielding a displacement 
. Thus, the creep rate, is by the following equations 


fe 
— 
| 


is 


CREE 
Here, m is the number of the dislocations located between the F- R source ant : 
_ the grain boundary, » is the distance br between two slip planes, are” 
defined by the identity in Eq. 1, k' ig the “specific rate of relaxation of an en- 
_ counter when f = 0, k"" equals 2 k’ sinh af, k is the Boltzmann constant, and T 
is the temperature Kelvin. The quantities, 1/a and 1/8, are called here the 
intrinsic shear stress and the intrinsic shear rate, respectively. F een 
- It is well known that the state of metals critically effects such onniuneeaa 
as cold working, heat treatment, annealing, etc. Such operations introduce fi 
certain types of bad sites into metals. Other types also exist. It is assumed 
herein that there are n kinds of flow units, each having a value of a and of f. 
Here a flow unit is characterized by the bad sites of an encounter such as aa’ 
Fo r the “ith ‘kind of flow: unit, Eq. 1 be written a as: 


"where, fj is the pa acting on the diffusing atom near — ‘ith kind of flow 


h unit. | Let xj be the fraction of the total area of the shear plane controlled xi 


it that each kind of encounter with same 


La rate, s, due to an appropriate concentration of stress on refractory encounter. ; 


The validity of Eq. 4 has been well established in polymeric systems by many — 
~authors.1 
- The pu purpose of the present paper is to apply Eq. 4 to the creep. of metals, i 


_ and to discuss several interesting problems related to relaxation theory. _ 


Bes Taikyue Ree and Henry Ind. Eng. Chem., Vol. 50, 1958, p. 1036. 
bytes 20 “Application of Ree-Eyring Generalized Flow Theory to Suspensions of Spherical a 
Particles. ae ae oon H. Maron and Perey E. Pierce, J. Coli. Sci., Vol. 11, 1956, 
i _ 21 “Application of nie -Eyring Generalized Flow Theory to Suspensions of renal 
- Particles. _ IL. Flow in Low Sear Region,” by Samuel H. Maron and Arthur W. — 
22 “Pipe-line Design for Non-Newtonian Fluids in Streamline Flow, ” by E. B, Chris- 
tiansen, N. W. Ryan and W. E. Stevens, Am, Inst. Chem. magne, —, 1, 1955, p. 544, 
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_ APPLICATIONS 


4: is first to creep p of phre and Al - Mg alloys. 3 
is found that the creep curves for these metals are well described by only one ‘= 
“ ~ kind of flow unit, that is by Eq. 1. In Fig. 3, the theoretical curves are eem~ 4 
4 pared i with the experimental cr reep data. 23 The values of 1/aand & are sum- 
marized in Table 1, with those obtained for other creep experiments. Table 1, 
- Col. 2 indicates the temperatures at which the creep (C) and the tensile (T) a 
tests were conducted. For the determination of the parameters, xj/aj and 
cs. occurring in Eqs. 1 and 4, reference may be made to the articles of Ree ae 
_ Eyring. 17,24 In Dorn’s paper,8, 23 ali the creep rates, §, are reduced as re- 


@ Tensile dote 422°K, 550°K, 650" 


activation heat, is 35.6 kcals for aluminum and its alloys. en gsi 


= "The application of Eq. 4 to the creep curves of the alloys,24 Al-Cu and Al- 
Zn, indicates that these cases cannot be expressed d by one kind of flow unit, but . 


require two kinds, that 


7 


a 


ing Structures in Alpha Solid Solutions,” by Oleg D. Sherby and John D. Dorn, J. Metals, a 
J 24 “Relaxation Theory of Transport Phenomena,” by Taikyue Ree and Henry Eyring, . 

Vet. Il by F R. Eirich), Academic Inc., New ‘York, 
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CREEP OF METALS 
+> sinh” 
— 


ar 


‘The th in Fig. 4 are calculated from Eq. 5 the 
ric values” mates in Table 1. By comparing the values of the 8; 8 (i= = 1 or 2). 


in their relaxation, (about five times as much as in the pure Al system), all 
_ the other is very easy, (requiring about ten thousand times less diffusion than — 


TABLE —CREEP OF ALUMINUM AND ITS ALLOYS 


Al-Mg (1 62%) 


422 477 528 
422 550,6 = 

©) 

= 
(Cc) 422° 
422,477 
422,550,650 
422, 


7) 422, 550,650 
422,550,650 


(3. 23%) 


| 


Al -Cu (0. 10%) 


422 
422 550,650 


422,477,528 


Al-Zn (1.62%). 


Annealed at it 650° for one hour. 


©) 422 2,477, 528 


408 5450 533 


-forthreehours. = 


~The creep data for the three aluminum samples listed ‘at the end of Table 1 


various grain sizes. 


is 


— 


ere 


Anneaiedst 10 1000 


were obtained by different authors: Al(2S-0), by Sinizer and MacGregor; Al 


- (99.995 %) and Al (2S-H18), by Servi and Grant. _ Oleg D. Sherby and Dorn23 7 


quoted these data to compare with their data. ‘Eq. 5 applies to the creep data 


of these aluminum samples with good results. ‘The parametric values ob- 


j 

Pe The following examples show the flow units are very sensitive to treat- — 

ments of the samples. . Dorn’ reported the experimental data obtained by de- r 


termining instantaneous creep r rates of the Al- -Mg alloys which were precrept 


toa given strain25 at 531° K. The authors apply ‘the following equation to the 


data: 


— 25 “Effect of Stress on the Creep Rates s of to gy Aluminum Alloys Under — 

Constant Structure,” 7 Oleg D. R. Frenkel, John E, Trans. 
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(0.55%) 
| 
| sa5 | | - =; 
— | — 
(28 - H18)° (C) 533 644.755 0.873 | 0.524— 2.70 | 0.02 
— 
— 
— 
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sinh"? $+— sinh 8. sinh” Bs 


ee parametric values listed in Table 2, are compared with the experiments. = 


& sees that the ; agreement is satisfactory. The sample of the aluminum base al- + 
= loy containing 3.1 atomic % magnesium were precrept to a strain, s=0.077 at 

_ 531°K under a given load, f = 10,000 psi. _ Other samples were also precrept & 
at 531°K, . but using different values of s and f; that ia, s= 0.049, f = 7500 psi 7 
Al-Mg (2.1 atomic %) and s=0. 077, f= 5500 psi for Al- Mg (1.1 atomic %). 

_ _‘The experimental data8 in Fig. 6 was obtained by using annealed Al-Mg al- 


and measuring instantaneous creep This case is described by 


WwW 
WwW 
w 
WwW 


Tensile dota 422°K, ,550°K 


° Creep Doto 422°K 
Creep Dot Dato = 
«Creep Dota 528°K 


+17, 800/T ($s IN I/HR. T IN °K 


4. —CREEP FOR THE ALLOYs: Zn AND Al-Cu. 


Eq. 6. By comparing the values for the present case (Fig. 8) with those for 
the precrept samples (Fig. 5), one sees that the annealed sample includes flow a 
“units which are harder to flow than in the precrept samples. : 
_ One may recall that the creep Curves for _ ordinary samples, which were a 
i obtained by measuring ‘steady creep rates, are described by | a single type of 
unit. This is marked contrast to the cases mentioned above, where three 
4 kinds of units were used. The authors believe that the difference is mainly — 
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M1 CREEP OF METALS a 
Dorn8 found from the experimental results. show in and 6 that that at 
~ high stresses, $ is proportional to exp (af) while at low stresses, it is propor- — 
tional to simple hyperbolic sine mile to oupress a creep 


13 a 


am 
Al -Mg(2.1 atm %) 


7 


_ dence of the -a@factor, which will be discussed later, was a riddle for many in- 

vestigators. these facts, Dorn® concluded that the simple 
; : - sine function proposed by Eyring16 is inapplicable for the Al and Al-Mg sys- | a 
tems. As mentioned the authors ’ theory involving more than 
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ne type of unit results ‘by 
Applying Mott’s theory of dislocation climbing, 26 Weertman!0 derived an 
equation for the bers of creep. _ He assumed that Lomer lockings act as ob- 7 


functions of Weertman’ s theory was applied to the creep of polyc rys- 


talline nickel by Weertman and P. Shahinian.12 The experimental points and 
7 Weertman’s theoretical curve (for 149° C) in their paper are reproduced in — 


1. full 7 is the -author’s theoretical curve which was | 


%) 


= 12.4 SINK! 4.88- 1075 


9 SIN '0.975° 


Al- -Mg (3. atm 


ealeulated, from Eq. 5 by using the parametric values: xj/ay = 4.98 x x 
= 8.25 x 10-10 (min), = 2.68 x 107 (dynes/cm2), and 
= 1.03 x 10- (min) Also in this case, the authors’ theory shows excellent 


. (Weertman’s 
= than does his theoretical curve shown here. ) Weertmant0 applied his a 


Mott, , Proc. Phys, London, 


“Dislocations, Plastic Flow and Creep,” N. 
Vol. A220, 1953, p. 


— 
q ( 
— 
ii 
— 
4 


"theory to the creep data of pure re aluminum (99.995%) obtained by — ~ 
Grant. — The agreement with experiment is not satisfactory, especially at high — 
stresses. _ The Servi- Grant data are very well described by Eq. 5 using = 
parametric values given in Table 1. The figure showing the agreement 


In a later paper,14 Weertman introduced the climbing mechanism due 
annihilation of positive | and negative leading dislocations on two neighboring “ 
& slip planes. The theory applies to the creep of hexagonal crystals, in which 
the Lomer locking does not occur. . His improved theory may apply to other 
a iy crystals also. - Although his new equation gives better agreement with experi- 


_ ‘ment at low : stresses, the deviations at high stresses are inevitable. La 


—+ 800°C - ~1100°C 


| 


:65,000cal/ mol. 


AD 


t(dynes/em*) 


7. CREEP FOR POLYCRYSTALLINE NICKEL, 


Intrinsic Shear Stress ond Shear Rate.- According _Dorn,® 8 the 

heat, AH and the intrinsic shear stress, 1/a, are independent of 

temperature. — “These two facts brought about the reduced creep curves as re- 
-_gards temperature shown in Figs. 3, 4, 6, and 7, The temperature indepen- 
dence of AH* is self evident, because at temperatures higher than 0.45Tm(Tm — 
| = melting temperature Kelvin AH* for creep equals AH* for self-diffusion,8, 27 

which is temperature independent. The temperature independence of 1/a is 

= known in high polymeric systems!7, 23,28 as it is in metals. ‘The ie q 


~ 27 “The Spectrum of Activation Energies for Creep,” by John E. Dorn, A Seminar a 
on Creep and Recovery, American Society of Metals, Cleveland, Ohio, 1957, p. 225. ~ J 
_ 28 “A Generalized Theory of Plasticity Involving the Virial Theorem,” Henry tell a 

and Taikyue Ree, Proc. Nat. Sci. U. S., Vol. 41, 1955, p. 118. 
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n Fig. 


ome the F-R source to the “The work done on each obstacle, 


Here, p (= o/r ) and are the surface density and the mole fraction of a 
sites, respectively, v (= a"), is the atomic volume of the creeping material, _ 
_ p is the probability of a site being a bad site yielding an encounter which © * 
total number, of these encounters of bad sites is calculated 


the ‘following considerations: The total number of bad site 


the counter plane is A p. _ The number of encounters of bad sites in the dislo- 
_ cation with bad sites on the counter plane is A p ¢, as the dislocation travels 7 - 


across the area, For a mathematical of this problem, see Ap- 

pendix 1. If the bad site, a’, vibrates energetically, that is the vibration energy 

than n hy(v= = vibration the avoids 


held 


a 
+ 


an 
Here, sate numerator is ‘the sum of vibrational states up to the nth level, and 
the "denominator is the total sum of vi states. Thus, one quinine Z 


_ Introducing Ec Eq. . 8 int into Eq. 7 yields: — 


Eq. 9 is equal to A As £/2. Thus ins: 


% 
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EM ‘CREEP OF META ALS 
If there are two kinds: of 
element , can ‘be represented by ¢ = + Here dp is mole 
fraction of bad sites in a pure metal, and ¢j is the mole fraction due to an al- 
ie _ loying element, while b is a constant relating the relative | efficiency of the bad 
_ site in encounters. If one is interested only in the range of Lon — accel 
a to Eq. 11a, 1/a@is independent of temperature, v being approximate— 
ly constant, and it increases with $1, the concentration of the alloying 
Both predictions accord withexperiment. 
From the definition of 1), one obtains the relation: 


ly constant. Thus, is independent of dy. 
The constancy of the values of Bin 1 Table 1, irrespective of varying are 


functional depeuinnce, a small change in a4 which “might arise with varying nde 

is practically w without effect on the results. ts. T Thus, the values of Bfor the Al- — 
ia, * Virial Theorem.—The temperature independence of 1/a may also be 

understood = ie 


“aie  - diffusing atom is assumed to execute a simple harmonic motion ina po- | 
pig tential trough. Let F and x be the force acting on the unit and the displace- 


ment, ‘respectively. to the virial theorem, the following equation 


_ ‘Here, KE is kinetic ‘energy, and the bars indicate the time average ¢ of the cor- 


responding quantities. The internal stress, f, (force per unit area) acts on the | 
atom, as the latter is displaced a distance, x, from its encom position. — 


Here G is the e shear modulus, and d is the distance between n two succ 


cessive 


} 29 “Introduction to Solid State Charles Kittel, John Wiley Sons, Inc. 


of alloying elements. This fact is interpreted by assuming that in Eq. 11b, Al 
 andA decrease with ¢,;. According to these assumptions, which seem to be 
— 
— 
| — 
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atomic planes. The quantity, a GA/27d, is critical stress, or yield 
stress, fo, .» at which plastic deformation occurs, a being a factor (<<1) which — ge 
= introduced to explain the fact that the observed yield stress, fc, deviates : 


stress, GA/2 md. Eqs. 12 and 13 yield 
Fx=aA 2%3 = k (14) 
ing to the calculations given in Appendix 2, the barred quantity in the 


second equality is 0. 1426. this result into the second equality 
Eq. 14, » one obtains 


In Table 3, Col. 2 the values of {, taken from references29, 30,31 are shown. — 

When there are no data of fc, the ultimate tensile strengths which always ex- 

‘TABLE 2.—CREEP OF PRESTRAINED Al-Mg ALLOYS 


i x x xi 
a q in. in. in hours: in hours hot hours 
| 1.00 | 0.500 

1060 — | 0.500 


a. atom % Mg)| 


1680" 
atom % Mg) 


fo, are used for calculating a from “The values of 


are shown in Table 3, Col. 3. The experimental values of ain Table 3, Col. 4 
are calculated fromthe slopes of the plots of Ins versusf, because Eq.1 is ap- 
proximated by * (1/28) exp f at high stresses. This method for obtaining 


units, Asa result, the a values obtained in this manner are for high stress. _ ; 


‘The ‘agreement of theoretical and experimental values of aare satisfactory. 


: Especially worthy of note is the agreement for the single crystal and for 'poly- 
a crystalline aluminum, although the f, values differ by a factor of 100. One 
= notes in Table 2 that whenever the ultimate tensile strengths are used the = ; 


; oretical values are smaller than the experimental ones. _ The po poor agreement 


_30 “Metal Reference Book,” by Collins J. Smith, Butterworth Scientific Publication, 
London, Interscience Publications, Inc., New York, 1955. — a 
my “Handbook of be and Chemistry,” Chemical Rubber Publis hing Co. » Cleve- 
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for gold is partly due to the use of the ultimate. for fe and 
partly di due to the fact that the enperimnatal value of a was calculated from a 


3.-CREEP OF PRESTRAINED Al-Mg ALLOYS 


Pt 109 2. 2.828 x 10° 
am 6.15¢ x 108 5.8 x 1079 9.210 x 105 
1.14 to 1.48% x x 108 77 x 10-7 


59) x 1077 


7.601 x 107 
Sn (Single crystal) 1.77% 10-7 


21> x 10-9 
4.459 x 10-3 


of Metals, Saul Ww. Dumbar and H. Huthsteiner, J. App. 


Earl R. Parker and M. W. Nathans , Technical Report, Institute of Ee cceanten Research 
7c’ _ Institute, University of California, Berkeley, California, Series No. 28,Issue No. 8, May > 
«High Speed Tension Test at Elevated Temperatures ,” by A. and Mg. 
-Manjoine, , Trans. Am. Soc. Mech. Engrs. Vol. 63,1941,p.A-77. 

“The Deformation of Gold Wire at Elevated Temperature,” by B H. Alexander, 
Dawson and H. P. King, J. App. Phys., Vol. 22, 1951, p. 439. 
bk “Creep and Recrystalization of Lead,” by A. A. Smith, Trans. Am. Inst. Min. 
ee. t «The Flow in Metals under Large Constant Stress ,” by E.N.da C.Andrade, Proc. 

AN “Creep of Lead and Lead Alloys,” by J. McKewn, J. Inst. Metals, Vol. 60, 1937, 
“Creep of Aluminum Single Crystals,” by J. Weertman, J. App. Phys. Vol. 27, . 
aa r “Some Fundamental Experiments on High Temperature Creep,” by John E. Dorn, — 
Symposium on Creep and Fracture of Metals at High Temperature, , Philosophical L Library, . 
_ New York, , 1957, p. 89; J. Mech. Phys. Solids, Vol.3,1954,p.85. = 
a —-™ «Creep of Tin Single Crystals,” by J. Weertman and J.E. Breen, J.of App. Phys. 
1 “Properties of Materials at High Temperatures. ” by J. Tapsel and W. J. 
Clenshaw, Great Britain Department of Science and industrial Research, Special Report, 
phen 'o “On the Variation with Time of Strain Rate of Metals under Static Stress ,” by A. 


_ Pomp and W. Lange, Mitteilungen aus dem Kaiser Wilhelm Institut fiir Eisenforschung zu 


curve which was obtained only over the range of 


= 
= 
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‘temperature range. 30, 31,32 ‘Thus, becomes independent of 
(Eq. 15) in agreement with experiment. Since is defined by the quantity, 
 AAQ 3/2 k T, the the independency is due to the linear increase of A Ag A3 with _ 
_ _ Hardness of Alloys.—One sees from Table 1 that the value of 1/a- for the 
aluminum- -magnesium alloys increases with the concentration of 
In Fig. 8, the values of 1/aare plotted against the concentration of | magnesium. 
bay straight line is the plot of Eq. 11, where v = 1.76 x 10-23, 6 = 281°, n=1, 
= 5.86 and d, = 8.95 x 10- -2. Over the range of low concentrations, the ex- 
aoe data fall on the straight line, while the point at 3.23% deviates ap-— 
preciably. Generally, the linear relationship (Eq. 11a) ) holds well for dilute — 
binary alloys, but deviates at high concentrations. 33, ial ef oo Sa 


The dotted | line in ‘Fig. is calculated from Eq. 1 the ‘ 


— 


A 


FIG. 8.-THE CURVE OF 1/et VS. CONCENTRATIO N OF MAGNESIUM 

 stresses”30 corresponding to the aluminum base alloys containing the atomic 

: fractions, d, of magnesium. | Here the ‘proof stresses are multiplied by an ar- 

4 bitrary factor of 1/3. 22 in order to convert the proof stresses into yield — 


& stresses. One sees that Eq. 15 fits the experimental data fairly satisfactorily. 


_As discussed above, 1/ais proportional to the yield stress. Thus, ultimate — 


Mia “Elevated-Temperature Properties s of Coppers rs and Copper r-Base Alloys,” by Clair 

Upthegrove and Henry L. Burghoff, American Society of Testing Materials, Philadelphia — 
J ie 33 «Physics of Metais and Alloys, ” by W. Boas, John ohn Wiley and Sons, Inc. New ‘York, 
atk, 


_ 34 “A Statistical Formulation for Creep of Metals,” by J. Lambert Bates, Taikyue 


and Henry Eyring, Technical 56, of Utah, Salt City, 
Utah, 1956. 
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METALS 


- stress, which are all functions of ald stress, increase with the concentration 
_ of additives in binary alloys. 30 Among these properties, hardness of alloys 
= been extensively studied. 35 According to the studies, hardness increases 
with additives in the low concentration range, passes through a maximum, and 
= decreases. 35 One / would expect the effect of additive concentrations on | 

— 1/a to parallel that for hardness. The authors, however, do not attempt here 


to . develop their theory « on 1/a over the whole range =e ene concentrations 


since the experimental data are not available. 


_ The authors wish to thank the National Science Foundation for their support &g 
of this work, The first authors expresses his appreciation to the Research 
. Committee of the University of Utah 1 for the the : fellowship granted to him. taal 
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Eq. 7 of encounters of two bad sites, a and a’ 


_and Fig. 2(a), is Ap. The number, A p¢, is obtained also from the follow- 
ww Atoms of two kinds, 1 and 2, of equal size form a close- -packed s space lat- 
tice , in which the number of nearest neighbors is r. Let Nj and Ng be the 

numbers of the respective atoms, the total number of the atoms being N (= Ny 


= Nz). T he atoms of the two kinds are assumed to distribute randomly, there _ 
¥ being no preference as between the paris, 1-1, 1-2, and 2-2. . The probability, — a 


W, fora a 2 particle of kind 1 forming the nid 1-1, in a given direction is: i) 


Dre (N, =1, i) C N,, r-i)/ C (Nel, r). 


‘Here C(N-1, r) = (N-1 + r) !/(N-1)! r! and indicates the number of possible 
"ways | of choosing r atoms out of N-1. Here the total number of available atoms © 
rine | N- -1, instead of N, since an atom of kind 1 is already known to be in a given 
; - position. © The product of the combinations in the numerator is the number of 
_ possible ways of filling the i positions of the nearest neighbors with atoms of 
_ kind 1, and the remaining r-i positions with atoms of kind 2. The probability 
for finding the 1-1 pair in a given direction is i/r for each complexion in the 
combinations, -1, i) C (N- -Nj, r-i). The ‘summation gives all possible 
_ ways for filling the nearest neighboring positions with atoms of kind 1 and 2. — j 


Using the N>? the 1 7) becomes 


_ 35 The Platinum Metals and Their Alloys,” by R. F. Vines, The International Nickel 
Company, Inc., New York, 
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Pica 
_ Thus one finds that the ee W, is @ iteelt as is, of course, necessary. — 
Since there are A, p particles of kind 1. on a plane of area A, the number of the — ; 
required 1-1 pairsisAp¢. Extending the above considerations, one canin- 
troduce into Eq. 1 18 the effects of the interactions between atoms and of the 


e the bar red quantity | in Eq. 14. T 


‘Here T and v are the pe riod and the linear | of the u 
-4 which is assumed to execute a simple harmonic motion. . There exist the fol- 
={@E-*) 
— Here, w is the Lee r velocity, m is the mass of the unit, g is the force —- #£#§#;j&§ 
. stant, and E is the total energy, that is, m xt/2. x= v 


| 
where d= Nj/N, that substi- 
a tution, i-1 = i, into the latter to: 
— — 
- 
— 
ie 
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19 into Eq. 20 


‘sin (me sin n 6) sin 6 a0 


"Here, the term, sin(7: sin n 6), is is in in the last 
being the Bessel function of the first kind of order 2k-1. . ‘The integration of — 
1) sin 6 dé from -77/2 to 7/2 vanishes for all values of k except k = 1, 
_ where the value 77/2. ° ‘Thus, Eq. 22 | becomes: 


{ad, (mh /2 = 0.1426 a . . 
where the value of 1 mi is obtained from a mathematical table. 
ie _ The above evaluation of I is an exact one. The following is an approximate 
> method. The quantity, v, in the last term of Eq. 19 is substituted by 2A/T. In- 
= the last term by parts, one obtains I = T= = 0.159. The deviation 


- of this value from that given by Eq. 23 is 11. 2h. as ae 


ween 


| 


— 
2) 
Let x=(a/2) sind. Then Eq. Stransformsto: 
— 
| 
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yee 1 M. ASCE and R. E. B. McCarthy, A. M. ASCE 


Experimental results on the behavior of three portal-type structural 


— 


steel rigid frames are reported. One such frame was subjected to oa 


al collapse of the frame. For me 


| 


variable loads applied correlation of the experimental lre- 


sults” with a refined was found to be. satisfactory 

ae of the normally employed procedures for determing the critical loads 

7 based on on the simple plastic theory appears to be conservative. ne se 


:" numerous engineering s structures magnitudes of the applied loads vary 
. _ independently of each other. Arbitrary application and removal of live loads, 
including wind and earthquake, is a common occurrence. Nevertheless, in the — 


_ _Note.~—Discussion open until June 1, 1960. To extend the closing date one month, a et 
written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the = 
_ American Society of Civil Engineers, Vol. 86, No. EM 1, January, 1960 - Seed 
a 1 Prof. of Civ. Eng., Univ. of California, Berkeley, Calif. 
2 Tank and Steel Co, Bakersfield, Calif. . Formerly, Graduate Stude 


— 
— 
— 
| 
nating plasticity inone of the members; and the third, to repeatedlater- 
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The of the limits of the various loads with 
arbitrary sequence of loading, termed “variable repeated loading”, as pointed - 
out by B. G. Neal3, was first recognized by Griining and Kazinczy. In an ex- 
cellent paper P.S. Symonds and Neal4 a systematic di discussion of ques - 


_ Two types of failures which may be caused by variable repeated loa 


s are 
recagiinell . In one type of repeated loading the reversal of stress is of such BP > 
‘magnitude that plastic flow occurs twice during each cycle. This situation 
a wy occurs when one attempts to break a wire by bending it back and 
_ forth. A possible breakdown of a structure at a critical section in this manner 
_ referred to as being due to “alternating plasticity.’ ” This type of failure is 
‘similar toa fatigue failure but : requires a much smaller number of cycles. 
‘The other kind of failure of a structure which is possible under a variable 
load is “incremental collapse.” | For sufficiently large magnitude of 
loads, after each cycle of load application an increment in the permanent or & 
inelastic deflectionof a structure takes place. Under the application of alarge 
number of loading cycles, large deflections develop. = 
= ‘For the simpler types of structures, theoretical studies of situations where 
alternating plasticity and/or incremental collapse based on ideal elastic- -plastic 
q behavior of a material occurs have been made.3 4 However, experimental 
q work to corroborate theoretical deductions has been very limited. 
Symonds? has found good agreement between theory and experiments on 


C. Massonnet y performed an extensive series of tests on two-span continuous 
steel | beams. _ The I- beams \ were NP n° 12 and the ‘spans each 114. cm 


“flections were re pointed out and good agreement with the theoretical predictions 


were indicated. At an earlier date K. Kléppel” in a ame series of tests 
reported no significant d deflections. 


A. -Gozum and G. Haaijer8 experimented with 4 WF 13 two - span con- 
tinuous beams, each span being 48 in. long. © Bending. took place around the 
weak axis of the member. To accentuate the effect, cyclic concentrated loads 
a applied at 0.4 of the span length from the center support. These expe 
_ ments indicated that beam deflection became stabilized under repeated appli- a 
- cation of loads which were somewhat higher than the calculated load at shake- 
(14.0 0 kips versus 13. kips). . At higher loads deflections continued to 


increase. 


_ 3B. G. Neal “The Plastic Methods of Structural Analysis,” - Wiley (1956), p. 269. 4 ~ 
a a B. G. Neal and P. S. Symonds “Recent Progress in the Plastic Methods of Struc- 
SP, S. Symonds “Comte. Loading on Small Frames” Final Report, 4th Congress IABSE 
y 


laminées”’ L’ *ossature Metallique (1954), No. 6, 318. * 


‘Kldppel “Beitrag zur Frage der Ausnutzbarkeit der Plastzitat bei 
spuchten Durchlauftragern” Final report, 2nd Congress, IABSE (1939), p. 77. 


8 ‘A. T. Gozum and G. Haaijer “Deflection ("Shakedown") of Continuous 
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evidently limited experimental work on steel structures subjected 


variable repeated load does e strong differences of opinion on this 
 matter.9, 10 On the other hand, 


= ESIGN N OF 
ne of ‘the of these experiments was to provide some in- 
‘ie on the behavior of simple portal frames designed for loadings en- _ 3 
countered in ‘practice. . Therefore, for the three frames selected for aS 


4 
FIG. 1.—LOADING CONDITIONS FOR FRAMES 
than the vertical ones. This is a condition encountered in the design of bents 
for wind or earthquake loading according to the standard building codes. Thus, 7 
_ in test No. 1 designed for proportional loading experiment, two vertical loads — 7 
_ were” set arbitrarily at4 W / each contrasted with a horizontal load of W as 
In test No. 2 the vertical load of 2 W was 1s specified to be applied in the | 
middle of the horizontal beam, Fig. 1(b). This load remains on the structure — 
to simulate the effect of dead weight as a horizontal force of ‘magnitude Wis be 
first applied inone direction andthen in the other. The four loading conditions 
constituting one cycle are shown in Fig. This experiment was to. 


Feat 


mut J - M. R. Horne “The Effect of Variable Repeated Loads in Building Structures De xs ’ 
signed by the Mastic Theory,” Proceedings IABSE (1954), 14, p. 53. 
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The frame for test No. 3, Fig. -1(c), was n made with unequal legs. The frame 
as thus proportioned inorder to achieve an incremental collapse failure when — 
_ Subjected to the repeated application and removal of a horizontal force Wz a a 
er Estimates of limit load for test No. 1 were made using the established pro- = 
_ cedures for simple plastic theory. 3 For the limiting shakedown conditions, 


that i is the largest the criterion fo: for r alternating 


or or neremental clase on 

m + Mmin p: 

_ denotes dies minimum moment computed o onthe elastic: basis, Myp i is the > mome 

at yield point stress, m represents the residual moment, and Mp is the plastic 

For the two cases of repeated loading, elastic analyses of 

‘ were performed and maximum and minimum bending” ‘moments for 
-Griie loading conditions were established. The fact that vertical loads: were 

continuously applied raised significantly the critical load for the case ofalter- 

‘The critical load for the incremental collapse case was cilia by apply- a” 
ing the principle of virtual work for an assumed mechanism. For this the — 
virtual work done by th the residual moments determined from Eqs. (2) must P 
All three frames were manufactured out of 4WF13 steel sections 
c . out, as shown in Fig. 2. Columns were pinned at the base. Corner joints mete ain 
a % made “rigid”. Note that the columns were cover-plated at these joints in order 


= to have plastic hinges form in the beams. sp avoided the necessity for con-— 


Frame No. 1, made as shown in Fig. 2, was loaded as shown in Fig. 1(a). 
The “proportional” loading was applied by means of three calibrated —— 
_ jacks. Pressure in each jack was independently controlled. A horizontal jack © 
4 was used to apply the desired horizontal force directly to the frame as may be 
- seen in Fig. 3. Vertical forces were applied by means of long vertical straps. in 
- _ Hydraulic jacks pushing against the main supporting frame and bottom yokes © 
attached to each pair of straps were e capable of deve loping the ne necessary ten- 
_Sile forces in the straps. (As may be seen in Fig. 3 a special adapter was — 
_ used under the right hand leg of the frame. This occurred only in test No. 3). 
an In order to prevent buckling of the horizontal beam during each experiment, 
an another channel section was placed in front of the frame as may be seen th 
‘Fig. 4. Frictional effects were minimized having narrow lugs attached to 
. _ A curve for the load W versus the maximum horizontal deflection of ll ie 
‘eee for test No. 1 is shown in Fig. 5. 
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for w versus ‘maximum vertical Deflection measure- 


the early stages of the experiment were made with ordinary | 
gages; later, with ‘machinist’ s scale. Horizontal deflections were measured 


located at four different positions on the frame gave inconclusive data i. 
Pity _ Note that during the later stages of -~ experiment that horizontaldeflection 


the ‘Shifting of the vertical loads, the column reaction on the right is ‘also ae 

7 creased. The severely overloaded column on the right, which bent and twist- 
ed somewhat, limited the carrying capacity ofthe frame. 

_ The second frame loaded as shown in Fig. 1(b) for studying alternating 

- plasticity yielded the load-deflection diagram shown in Fig. 6. - Note that for 

the first nine experimental points several cycles of loading were applied at 
each load. The number of loading « cycles is indicated on the diagram by a 


number at each experimental point. The last six experimental points are 


- single applications of the horizontal force W acting to the right, together wit 
ot The response of the frame under repeated cycles of load application may — 
. The average strain for gages at 1 and 3, located 


_ “shookdown” for these loading conditions. “This is n not the case for Ww equal to s 
6,200 lb, where strains somewhat increased with each. cycle. Certainly at W 
; —_— to 6,800 lb, the structure showed no evidence of shaking down. atta 
AS may be noted from Fig. 6, large horizontal deflections were observed - 
when large non- “cyclic loads were applied to the frame. . The final carrying 4 
aa of the frame, as in the case of test No. 1, was limited by the capacity 
s the right-hand column. Axial load together with bending moment, both mag- 4 
Bg nitudes augmented by horizontal deflection, caused severe bending and twist- 
in the column. ‘This may be seen by examining Fig. 8. 
Thethird frame was tested to furnish information on incremental collapse. 
The load-deflection diagram for this frame is plotted in Fig. 9. For each of . 
_ the first twelve experimental points numerous loading cycles were applied. 
i Fig. . 9, numbers by the experimental points indicate the number of 8 3 
at that load. Each of the last two experimental points was obtained for one oo" 
load application with all forces acting onthe frame. 


aan = Fig. 10 strain reading at gage point 1 for various conditions of loading 
1) number of cycles are plotted. From this diagram it may be —aan ; 
_ that for W up to and including 6,800 lb the frame shakes down. At higher loads | 
_ the strain appears to become progressively larger with each additional cycle. 
_ The same conclusions on the behavior of the frame may be reached by examin- 
. ing the horizontal deflection-cycle diagram plotted in Fig. 11. Excellent corre- 
lation between the two diagrams of exhibiting shakedown trends is evident. 
os After the experiments with the frames were completed, suitable sections — 
___ Were cut out from the frames. - These were used to check the dimensions of 


= rolled members. No serious discrepancies with standard dimensions were 


it was found ‘strength of steel should be at 37. 5 ksi. 
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“a The principal results of these experiments are summarized i in Table ta 
ies eperene the experimental loads We some vagueness regarding their actua 


=e 


—STRAIN- CYCLE DIAGRAM FOR ALTERNATING | PLASTICITY 


i values exist in case. The experimental “ collapse load” for the propor - 


tional loading test actually has been based on a deflection criterion. 


taken to be ten times the normally criterion for 
if 
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working g loads of 1/360- th - the span length. - From Fig. 5 it should be noted 


_ that the horizontal deflection was used as a criterion for d deflection. . If i maxi~ 


In test No. 2, the experimental load We for alternating plasticity was set at — 

- 00 kips. As seen from Fig. 7, this is conservative. It is conceivable that at 
ro loads as 6,200 lb, 6,400 lb, or, perhaps, at 6,600 lb the structure still 
would shake down. Howev er, experimental data are lacking to make this con- | fae 


clusion. Moreover, “since for alternating plasticity there is. some reason 


4 
| 


> 
4, 


ee which shakes down, the choice of 6,000 lb load appears to be p proper 
if basedon the curve of Fig.6. From Fig.7 note the erratic behavior of strains 
at W = 6,800 kips. . Also note that at total of 75 cycles of loading was applied to ~ 
- the fi frame without causing failure. In the last stages of the experiment the 
Sy 
frame performed satisfactorily under the application of proportional loading 


. pen In test No. 3 the incremental collapse load Ww was taken to be 6.8 kips. 
2 this load, as seen from Fig. 10, lita at gage No. 1 shook down. ‘This was + 
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a at least a as evidenced for the number of cycles performed i “. these aaa, 

_ Again it should be noted that after applying 204 cycles of loading, the frame 
performed well under the application of large proportional loadings causing 
py _ large amounts of deflection, Fig. 9. Note that the linear relationship heen 
7? load and deflection ceases to exist somewhat before the critical incremental 


c 
1, results of the conventional analyses are 


ae In the fourth column of Table 
epee. _ The reported values are based on calculations for frames with the © 


centerline dimensions and loadings as shown in Fig. 1. The yield ‘strength “a 


i used in calculations was 37.5 ksi and standard dimensions of cross sections 


TEST “3 
12.—ASSUMED COLLAPSE MECHANISMS 


Calculated values of W based on conventional analysis: are low in 
- son with the experimental results for all three frames. | The amount of dis- 


"agreement on a percentage basis is shown in Table 1 , Col. 6. “In general, , the 
amount of disagreement on this basis is substantial. - An improvement in the rs 
agreement can be achieved by making a more refined analysis of the frames. i. 
In this refined analysis, the Kinematic collapse mechanisms instead of being © 
- taken as shown in Figs. 12a, b, andc, (conventional analysis) were re assumed _ 
to be as shown in Figs. 12d, e, andf. This is based on the assumption | that — 
. plastic hinges form in the beam at the inner face of the column to beam con- 
Be nection. As may be Seen by examining Cols. 6 and 7 of Table 1, , aconsiderable = 
_ improvement in the agreement of theoretical and experimental Tesults is thus _ 
a obtained. This shows the sensitivity of a frame to a minor | change in the lo- q 
cation of the plastic hinges on the load. 


Semee = reproducible under repeated cycles of loading as may be seen from Fig. 1. i 
| 
EST 
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FRAME STABILITY 


Based o on these several conclusions be reached. 


The most important ones are: 


(a) All the experimentally determined loads were higher sidiaidietaiaain ie. 
Strain hardening of the material is undoubtedly y an important ‘contributing - 
= in causing this disagreement. In situations where numerous loading =. 
cycles are applied, the effect of strain hardening may be more pronounced than ia, 
_ that encountered in one application of aproportional load. However, the | agree- » 
- ment achieved between the refined analysis < of frames and experiment may be | 


"considered reasonably satisfactory. 


— (b) The importance of considering deflections of frames for determining 


their true carrying was strongly brought out these 


“(type i  Conven- 


| 


‘Proportional 
Variable | 
Repeated 
(alternating 
plasticity) 


Repeated | 
(incremental 


collapse) 


a Limit loads for ‘proportional loading based on the conventional analysis for tests 
numbers 2 and 3 are 6. 54 kips and 6.18 kips, respectively. the refined analy- 


sis, they are 6. 4 Kips and | 6. 34 34 kips, respectively. 


This contrasts with the situation normally foundin the experimentswithbeams. _ 
The large horizontal deflections of the frames tested caused “secondary” bend- aa 
_ ing moments in columns amounting to some 50% of the primary ones. of In de- 
_ termining true ultimate ‘capacities: of frames such effects must b be - duly con- ae 
Excellent performance of steel frames under loading after 
fs considerable number of cycles of repeated load applications should be noted. _ 
Further work in in this direction i is needed. , Experiments | for determining failure 
when alternating plasticity occurs would be desirable. Information on the ap- aa 
plication of a large number of loading cycles to frames in situations where 


incrementalcollapse occurs also would be of interest. The question of the q 


_ fects of previous load applications on the subsequent behavior of structures 
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The experiments in n this. paper as pilot tests. 
No report of experimental work with welded rolled steel frames under the — 
of variable repeated loading is known to the authors. However, such 
information is needed for correct plastic design of srdateid frames for wind 

and earthquake loadings. The results to date are encouraging by indicating 7 

conservative nature of the conventional methods of plastic analysis. 


: mt The silts wish to acknowledge with gratitude University of California ar and 
National Science Foundation (Grant NSF - G- 5953) for the funds made available 


ic 
* conducting and reporting this research. 
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LOADS ON RIGID-PLASTIC BEAMS* 


By P.S. Symonds! 
ink 


= he problem is considered, ofa or i beam crossed is a load 


is too large for the beam to support statically. The ‘main question 


fs considered is enettine or not there is a critical magnitude of the travel- 

7 


load, shove: which the crossing cannot be made at any finite speed 


causing total plastic coll collapse. a critical “magnitude: was 

predicted in an analysis made by E. Ww. Parkes. The present study in- op 


dicates: that a critical | load will not be obtained if the mass (of the the beam 


E. W. Parkes has recently the problem a load 1 
‘ideally plastic beam, 1 load being greater than the load under which 

the beam would undergo plastic collapse if it were stationary at the beam 


= the midpoint, the deflection increased irrespective of 
 Note.—Discussion open until June 1, 1960. extend the closing date one month, 
= ‘wrtaen request must be filed with the Executive Secretary, . ASCE. This paper is part x 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the 
Society of Civil Engineers, Vol. 86, No. EM1,January,1960. | 
_ 8 The work reported in this paper was partly supported by Contract Nonr- -562(10)be- “4 
_ tween Brown University and the Office of Naval Research, Department of the Navy. __ 
= Prof. of Eng., and Chmn. Div. of Eng., Brown University, Providence, R. I. — or 
2 prof. of Civ. Eng. University College of Swansea, Wales. Visiting Prof. of Eng., 
W. Parkes, “How to an Unsafe vol. Nov. 
pp. 606- : 
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stant speed, Parkes’ analysis showed that as the travelling load approaches 
a 


He therefore that the travelling not ‘exceed 


this regardless of its speed of crossing. 
‘This result is surprising. With the assumptions ‘made, one would expect 
that a load of any magnitude could cross the beam if it travelled fast enough. : 
% _No physical explanation is apparent for the existence of a ‘critical- load ma : 
assumption that the beam is massless seems at first a reason-— 
SS one for the case in which the travelling load is large compared with the 
beam weight. a However, difficulty arises with a massless rigid-plastic b beam 
— to the fact that no kinetic energy can be stored in it. it therefore cannot 
be in motion after the travelling load has ‘crossed, since this would involve 
3 the expenditure of energy in plastic deformation which, in turn, must be <i 
q anced by a decrease in kinetic energy. ‘Nevertheless, for the case considered, 
Ee one would expect that if the travelling load exceeded a certain magnitude the 
a 7 kinetic energy created during the passage of the load could not all be ab- | 
m 86§6 sorbed by plastic work while the load was still on the beam, so that above this j 


magnitude further deformation would take place after the load had gone off 

the beam. It may be expected therefore that when a massless beam is treated, 

_ difficulty will arise unless the moving load has a magnitude such that all de- 

_ formation is completed before the load reaches the end of the beam. _ This is 
evidently | the source of Parkes’ result: that deflections s approach infinite mag mag- 

¥ nitude as the travelling load approaches acritical value. 

EA In this paper we write the equations that would have to be solved ina gen- 
eral treatment of Parkes’ problem for a rigid-plastic beam with mass. We 
extract some information from these equations, but do not solve them. In-_ 
stead, we give thesolution of a related problem, namely, that of a rigid- -plastic 
beam strengthened in such a Way | that plastic | deformation occurs only at the 

Soa It is found for the problem | treated that there is no critical a 
when the beam has no mass. - This example provides at least a strong indica- 
_ tion that for any beam with finite mass there is no limit to the magnitude of _ 
the travelling load; for any given magnitude of the load the deflection can be 
limited to any pre-assigned size, if the speed of crossing can be increased to 


ti _ The analysis is again of the rigid, ideally plastic type, so that curvature» 
pe change can occur at any section of the beam only if the bending moment there 
has the magnitude Mp, which is taken as constant and equal to the “limit mo-_ 
ment” or “ fully plastic “moment” of the plastic theory of structures. This a 
of treatment may be accurate when the plastic deformations are ‘sufficiently 


5, 6 compared to elastic ones. We ‘shall not attempt here to take account 


Suppose the loads on the beam are as shown , consisting of 
“EST, 


static uniformly distributed load W= p pt g and a travelling load P= 
where p is the mass per unit length of the beam and m is the mass of the — 


‘ _ moving weight. Under W alone the beam would fail by plastic collapse at a_ 


_ 4£E.H. Lee and P. S. Symonds, “Large Plastic Deformations of Beams Under Trans- 
verse Journal of ‘Applied 1 Mechanics, Vol. 19, 1952, pp. 308- S14. 


a P.S. Symonds and G. F. Leth, “Impact of Finite Beams of Ductile Metal,” » Journal 
of Mechanics and Physics of Solids, Vol. 2, 1954, pp. 92-102. 


Bit E. W. Parkes, “The Permanent Deformation of a Cantilever Struck ee er 


at its Tip,” Proc Royal Society, London, Vel. 238A, » 1955, 5, PP. 462- 
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load 8 If in addition Ww the | force P is applied ai at the midpoint, 


collapse will occur when P has them magnitude % ‘such that > Po Mp. 

In our a ia it is convenient to define the loads by | means ofr ratios a, B as as 


is the magnitude of P pont that plastic collapse occurs when both. loads ee 

applied together, P being stationary at the midpoint. The travelling a 
7 _ located by its distance s =v t from the left-hand support, t being time and 

= 8 the velocity, here assumed — derivatives sre 


= j 
Motion “rigid treatment), when at some 
cross section the maximum moment in the beam reaches the magnitude Mp; 
and plastic-hinge action can then take place at sucha section. In the subse- 
quent motion there is a plastic hinge at a cross section x from the left-hand | 
- support, 4 The segments to left and to ) right of the hinge move as rigid bodies 
" instantaneously of length x and & - x, respectively, and with angular velocities — J 
9 and ¢, respectively. As the deformation proceeds the hinge moves, in gen- e 
: eral, leaving behind it a a zone where permanent curvature has been acquired. 


at the left and at the right-hand supports, respectively . The transverse ve- 
locity of the beam is continuous at the hinge section, so that the angular ve- 7 
locities 8 and @ are related by (£ - x) d =x 6. When the hinge is moving, 7 


ever, the is not continuous, and the following relation 


; _ In Fig. 2 are shown free e-body diagrams of | the beam, including the given 
loads P and W and the d@’Alembert forces corresponding to the travelling mass _ 


m and the distributed mass p per unit length of the beam. The tr avelling ; 


— 
— 
— 
— 
— 
= 


2000. 


is to stick to the beam, so that if it is to the left of of the 


4 


_ Fig. 2(a), the net effective force is P - ma-= P-m (so +20 s 6). IfP is to 
the right of the hinge, Fig. 2(b), the net force is -m 

2s¢]. We assume throughout that displacements are so small that dis-— 
tances along the beam can be replaced by their horizontal projections). pal ‘= 
_ Equations of motion are easily derived from the free body diagrams a 


Fig. 2. The forms given below were obtained by the method of virtual work x 


pair was derived using virtual displacements of ‘the type indicated by Fig. ¥ 
while the second in each pair used virtual displacements as indicated in Fig. * 
8b. In the system of Fig. 3a the only internal work is done by the bending — 
- moment Mp at the hinge, while in the | 
be by the shear force at the hinge | section. "This shear 


p x36 


subject to starting conditions for each interval. 


_ Deformation first begins when the load reaches a position Sp anda plastic. — 


hinge forms at. These obtained by sctting the 


: hand sides of Eqs. 3 equal to zero, as follows: 
So 


— 
co n on the bending mc 
* (at ls, P=me=0 
— 


Here and the non- -dimensional co-ordinates = and o = s/t 
will be used to locate the plastic hinge and the travelling load, respectively. oar 
Since o = v it is convenient to use this as non- time aca scale 


‘By.a limiting process, writing = 09+ € and assuming that § = 
“valid for small € € the initial lal velocity o of the hinge point found from Eqs. 


— — 
— 


Motion when the relative angular at the zero, 
hence when @ = 0 in our case. This can occur either while the load P is 
still on the beam or (for sufficiently large loads) after P has moved off the 


where =v is used as the independent variable and differentiation 


At the instant when motion hand ‘side of Eq. 
corresponding final e of is found to be 2° independent the | previ- 
ous history of deformation. If the load magnitudes are such that t the deforma- 


“hinge | is re 


Since we are interested in the question of whether or not. there isa limiting 


_ value of the ratio P/P, _ rather than in the d details of a complete solution, we 


not integrated general equations. Instead, an ‘approximate solution be 


found by putting plastic hinge arbitrarily beam midpoint, 


the motion. Thus Eqs. 3a, 4a, and 5a were solved with = 


errors by this are of magnitude , but it 


g seems likely that our approximate solution will correctly show the qualitative 
aspects of the complete solution ‘relating to the question of a critical magni- 
a tude of the travelling load. It can, 1 of course, be interpreted as the correct SO- mio 


we have = = , and 3a, , and 5a can be as fol- 


P/w = =a) a) 


AS — 
— 
— 
aa 
— — 
—. 
ad 


4 ‘LOADS ON BEAMS 


» 


ti are @=6=0 when = Oy, where % a obtained by setting the right- hand 


of Eq. 10a to zero. Since simple result is obtained 
that 


-1 


122 12m _ 6 (i - - 
6=0 wheno =o Of where is found by setting the right of 
Eq. 12a equal to zero, the final angle @¢ is found by putting o = a in Eq. 12b. : 
On the other hand, if @ does not reduce to zero in this interval, motion con- 


tinues after Pp ‘has left the e beam. The final angle Og is then given by : 


me 


it! is found by examining Eq. 12a that if 1 .293 the beam 
_ continues in motion after P leaves the beam, since then 6 (1) > 0. By Eq. 4 & 
this means that p = P/P, must be larger than 1.707 in order for motion to 
persist after o = 1 , irrespective of the mass of the beam. 
We can now look» into the question of whether or not there isa finite limit 
—— the load P that can cross the beam. hs Examination of Eq. 14 shows that for f 


— 0 there is no ) such limit. For very large v values of P/Pc we have % = 


—- 
Ace 
a — 
— 


2 
PB 
/6B (1 - a) in terms of the basic ratios a 
and Pt/4Mp tt it is that if a > 0 quantity 
g becomes very large only | when P/P. becomes very large. Evidently 
| "therefore, if the beam has mass , an arbitrarily large load P can cross the A 


beam without causing an excessive deflection, provided its crossing speed V =v 
Fig. 4 are shown curves of v2 asa of P/P, for the values 
As an example, if a=7, @£ = 50 ft, and the maximum deformation — 
of is to be held to 0. ‘ an overload factor P/P, of ae can be applied provide = 
crossing is made with a speed of about 81 fps. 


markable result that the and deflection become infinite as the 


‘overload B= =P/P¢ approaches 1 1.707 In the | ute 1, 


| 


The final deflection is i 


of Eq. 17b shows that as" go 71 - 0. 293 from above, (that 
‘1. 707 from below) the » value of of “approaches unity. As 
happens, the. equations | show that the maximum angular velocity and the 
final angular deflection 6¢ both approach infinity. Evidently therefore the load» a 
is ratio P/P, =1 107 is a critical one; our theory predicts, , in a manner similar 
that arkes, that the vehicle a weight equal to or exceeding 1. 707 
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"LOADS on ON | BEAMS 


plastic beam with weak mid- section. - The fact that a critical ‘limiting load is : 
_ found when the beam is assumed to be massless, while no such critical load 
found when the beam has mass, is the significant result. 

iat we had considered a uniform ‘beam the solution would have’ involved a 


to zero. Thus in the case of the massless beam treated by Parkes the hinge © 
ie If the beam is supposed to have e any finite mass (however small) the hinge 
will lag behind the load by a finite distance . Then there is no difficulty about | 
7 continuing the plastic deformationafter the load has left the beam. As soon = 
load reaches the vight- Support the — moves back ‘the 


the mass of the beam decreases to zero wanes amount of the lag also | decreases 


— 
will produce in GELLCCUONS OF the if it crosses with any finite 
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center, where! it is when the deformation comes to an end. 
: the other hand, if it is assumed at the start that the beam is massless, then | 
the hinge must cross the beam at the load point. This leads to no difficulty so 
long as the deformation is completed before the load reaches the end of the Ey 
a beam. But since there is no mechanism by which the deformation can continue 
after the load has crossed the beam, the load magnitude at which deformation 
os just completed as the load reaches the support must appear as a me 
load value that cannot be exceeded. Such a critical load is evidently a a special — 
‘property ofa massless rigid- plastic beam. 
a both this treatment and that of Parkes, the basic equations were derived — : 
on the assumption of very small deflections , 80 the prediction of * infinite” de- 
flections is not to be taken literally. Nevertheless a small- -deflection theory 
a =, well predict a physical significant critical load (as for example in Euler ro 
buckling theory) and this is the question under examination. If it was not as- A 
sumed, in deriving the differential equations, that angular deflections and 


“velocities were very small, then additional terms — appear in the ee 


cle, ‘development of axial forces, ete. In most cases the equations would prob- _— 
F ably be difficult to solve « except by ‘numerical methods. Moreover if the 4 
was massless the fact would still remain that kinetic energy could not be 
stored in the beam after the load had crossed, It is therefore probable that a 
critical load would still be found, corresponding to the load value below which oa 


the deformation would be completed during the time the load was crossing the 
An analytical anaes. can in fact be easily found for one case, with ‘results| : 
valid for arbitrarily large and 6; this is the case in which v_ stands for the 
(constant) horizontal component of velocity, the beam is held by fixed pins at 
_ the supports, and no axial forces are assumed to arise because of large de- a 
“flections. It can then be shown that Eqs. 16 give the solution for the massless © > 4 
beam witha plastic. hinge at its n midpoint, provided that in Eq. 16a 6 is i 4 
placed by 6/cos2 @ , and in Eq. 16b @ is replaced by tan @; in both equations 
i? = Stands for the horizontal ‘component o of displacement from the 


an eccentric axial compressive loa load which is at a value | P for a finite 
“ _ interval of time. If the strut has mz mass s the load P can reach any size; ; provided 
‘the application time has a sufficiently small maguiude the maximum deflec- 
: tion can be limited to any given small amount. - _ However, in cases 3 where the 


7 Jules Verne, “Round the World in Eighty Days,” Chapter 2 28. 
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Vided large enoup peed OF Crossiig Could be attained. Ihe present 4 
a ac analysis tends to substantiate this view as far as girder bridges of ductile = 


‘metal are concerned: failure by the mo: most form 
of failure in this case, and the rigid-plastic beam analysis then contains the — 

essential features of the problem. - The analysis brings out the point that the 


mass of the beam must be taken into account, as otherwise an upper limit is | 
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RING MECHANICS 


_F roceedings of ‘the American Society of Civil Engineers 


Load carrying capacities 0 of plates ‘subjected to axially ly 


trate the use of the charts” the ots a vessel 


“ ie — plate is an important structure which is is widely used in ma- 


chine constructions, | ‘such as engine pistons, load speakers, heads s of pressure 


vessels, etc. In the} design of the circular plates the allowable load which has 

_ been conventionally determined on the basis of elastic analysis is governed by 

of maximum deflection and/or maximum stress. the 
£ case in which the ‘strength of circular plate is of primary concern, 
elastic analysis” gives ‘the maximum allowable load as the maximum load at 

_ which the stress state at no point violates the yield condition Specified by a 


"highly 4 desirable to raeeton a circular plate taking the plasticity of the material 3 


‘into consideration to achieve a better utilization of material, and to obtain a 4 
lower strength-weight-load ratio. In general, an exact plastic analysis of | 7 ; 
‘structures would be a formidable task for designers. However, the method of a 


a. Note.—Discussion open until June 1, 1960. Separate Discussions should be submitted 
Se _ for the individual papers in this symposium. To extend the closing date one month, a : 
_ written request must be filed with the Executive Secretary, ASCE. This paper is part os 
_of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the ra 
Society of Civil Engineers, Vol. 86, No. EM 1, January, 1960. 
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Selected theory Of failure. in other words, under such a load, the plate a 
ES saa elastic. However, a circular plate is capable of sustaining a load higher tar i 
4 the maximum elastic allowable load if a nermissible plastic deflection rather 
— 
— 


by R.M. Haythornthwaite.? The effect of the membrane stress 


Haythornthwaite and good agreement between the analytical results and test 
data was found.8 In 1955, Hopkins and A. J. Wang reported their studies on the 


Limit | design provides a simple yet reasonably satisfactory approximate method 
* neglecting the strain-hardening of the materials. 2 Hence, design of circu-— 

plates based upon the plastic limit load has been suggested. 


_ The elastic analysis of a circular plate with axially symmetrical loading 


has long been known and is given intextbooks on elasticity and strength of ma- as 
application of the elastic analysis to circular plates with various 


types nf loading and boundary conditions was made by A. M. Wahl and G. Lobo 
in 1930.3 Only in recent years, however, has the plastic analysis of circular _ 
plates made great progress. Attention has been ‘essentially confined to the de-— 
termination of the load carrying capacities of circular plates subjected 
axially symmetrical loading. The first eyo work on this subject was due 
to H. G. Hopkins and W. Prager in 1953,4 _ although the concept of upper bound | a 
collapse load of plates was developed by K. W. Johansen9 as early as 1932. 


% cular plate by employing the concept of Tresca’s yield condition and the asso-— 
ciated flow rules. In their analysis, the circular plate was assumed to be iso-— 
tropic and rigid-perfectly plastic. _ The case of a circular plate with ith overhangs | 
was investigated by D.C. Drucker and H. G. Hopkins.6 The load carrying ca 
pacity of a circular plate with a rigid center boss was ; studied analytically and 


on the load carrying capacity of circular plates was analyzed by E. T. Onat and 


_ load carrying capacities of circular plates of perfectly plastic material with — 


* arbitrary yield conditions. ail The determination of the load carrying capacities 


_ of circular plates of anisotropic materials was given by A. Sawczuk!0 and Le 


at, 
“Plastic Design Methods —Advantages ‘and Limitations ,” Proc. 


Society of Naval Archite ects and Marine Engineers , pp. 172- -196, 


Wahl, M., and G. Lobo, “Stresses and Deflections in Flat Cizoular with 


- Central Holes, > Trans. ASME, Vol. 52 (1), APM- -52-3, Pp. 29-43, 1930. pe 


Hopkins ,H.G. and W. Prager, “The Lead Carrying Capacities of Circular Pl 


Appl. Mech. and Phys. of Solids, Vol. 2, 
Johansen, K.W., , “Moments of Rupture in Cross Reinforced Slabs,” Proceedings of 
International Association of Bridge and Structural Engineering, Vol. 1 »Pp. 296, 1932. 


6 Drucker, D. H. G. “Combined Distributed Load 


bie 33 on Ideally Plastic Circular Plates ,” Proc. 2nd U.S. National ie cnn of Applied Me- an 


Haythornthwaite, R. M., “The Deflection of Plates in the Elastic-Plastic Range,” 
Pre: 2nd U.S. National 1 Congress of. Applied | Mechanics, PP. 521- -526, 1955.00 a 


8 Onat, E. T. and R.M. Haythornthwaite, “The Load | Carrying Capacity of 
Plates at Large Jour. ar. Appl. Mechs., » Vol. 23, 55, 1956. Dingo: 
3 
or Hopkins = ¥ G. and. A. J. Wang, , “Load Carrying Capacities | for Circular Plates of ‘yn 
Plastic Material with Arbitrary Conditions ,” Jour. and Phys. of 


10 Sawezuk, A., “Some Problems of ers Carrying Capacities of Orthotropic and a 


| Non- Homogeneous Plates,” Arch, Mech, Stos., Vol. 8, No. 4, pp. 549-563, 1956. ae i 


11 Hu, L. W., “Modified Tresca’s Yield Condition and Associated Flow Rules | for 


- Anisotropic Materials and Applications,” Jour., Franklin Institute, Vol. 265, pp. 187- 
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-CIRCUL LAR PLATES 
In this paper, the —_ limit loads of circular plates for various on of 


and boundary conditions are determined. Accordingly, design charts 
for various cases are constructed. 


A 


plane and points downward. “Let r dentin the distances of points in the middle 

plane measured from the center of the plate, w 
principal curvatures K and Ko are: 


The radial bending moment and the ta tangential moment are 
to be positive whenthey produce tensile stress in the bottom fiber of the plate. 
s All the forces are considered to be positive when they act downward. _ From <4 
the consideration of equilibrium of an element of the plate, the following equa- 
tions of equilibrium are obtained: 


the My and Mg and the shearing force are 
: i intensity per unit length, and the external load p is defined as intensity per 7 
unit area. The lower limit of the integral in Eq. 2bis eee ce by | the bound- 
ary condition of the problem to be considered. 
= When the circular plate is elastic, ‘the bending moments are are ) related tw am 


rincipal curvatures the manner: 


Where D= 


— — 
— wi” 
—- 
loading. The undeflected middle plane is in a horizontal position. Passing 
| 
| 
i — 
— 


in in elastic 1 range can The ‘procedure : and the solutions to 
various cases can be found in textbooks on elasticity and strength of materi- - 


a plastic region will form in the circular | plate e when the stress state in the 
== the critical state specified by. yield condition. In this paper, 


Where Or radial % is the tangential stress” and is the yield 


_ stress is tension of the material under consideration. The plastic flow will be - 
a governed by the associated flow rules which state that the three principal plas- - 
tic strain rates at a point are proportional tothe three direction cosines of the 
! . exterior normal to the yield surface at the point representing the stress state 
_ of the material point under consideration. — For a point on the edge formed by — 
-, two. adjacent faces, each of which has a unique normal at this point, the princi- 
oa pal plastic strain rates are assumed to be proportional to alinear combination | 
of the two sets of direction cosines of the two normals. = 
ae For the determination of the plastic limit load of circular plates in which > 
; the bending stresses are of major concern, itis assumed thata perfectly plas-— 
- tic material will deform plastically only after the state of bending moments 
reaches a certain critical combination as specified by the yield condition. The 
Limit 1 moment M, which will cause a beam of unit width to become completely — 7 
BF aaron is taken as a basis for formulating the yield condition aforementioned. © 
Consequently, for a circular plate of rigid-plastic material, deflection becomes 
possible only after the bending moments ts throughout the whole plate o ora part of of 
reach the yield condition defined by: 


Where Mo = and is the stress in tension. The principal rates of 
deformation are governed by the associated flow rules as as vel’ in in Fig. 1 for 
determine the plastic limit load in circular plates, the  distribu- 
ten in elastic state must first be known. | Then it is necessary to investigate e 
propagation of the plastic zone. The moment distribution in the plastic zone 
must satisfy the equations of equilibrium and the yield condition. Since the 
| plate is considered to be rigid-plastic, the incipience of deflection is governed © 
sia by the associated flow rules and the boundary conditions. The plastic limit — 
load is reached when a part of or the whole plate ‘becomes plastic and the in- 
_ cipience of deformation is found to be possible. Then the idealized circular — 
plate will collapse under the limit load so et 


BENDING OF CIRCULAR PLATE BY SHEARING FORCES 4 
To the general procedure of determining the plastic 1 limit load 
of a circular — let us consider the case of bending of a circular plate > me 


— 
3 
— 
— 
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P First the moment distribution in the plate inthe elastic c range can be deter- z 3 


(6b) 


7 es these results, it is checrved that both the tangential and radia 


(Mo, Mo) 


chy. 


= 


Mr= = 
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moments are positive throughout the ; plate. . Furthermore, the tangential mo- 
ment at the inner boundary is the greatest. It can be concluded, then, that the ; - 
inner boundary first becomes plastic when the shearing force Qh reaches the | 
= -— eritical value. | a That is, the tangential moment at the inner boundary reaches © 
the full plastic condition first. The stress state at the inner boundary is rep- 
ey resented by the point Cc of the yield condition, -as shown in Fig. 1 ‘since t there 
As the intensity of the shearing force ce increases, the plastic region will © 
r and from the inner boundary outward. However, it is necessary to know 
whether some new plastic region may he formed during the expansion of the 
plastic region already developed, because the load carrying capacity is 
mined by the of deformation which in turn is 


laste Region 


a us consider that the plastic region has extended from the'tener 
 ealec to radius c. Namely, the portion of the plate within the circle of | 
_ radius c is” plastic, while the portion of the plate outside the circle of radius 
- ¢ is elastic as shown in Fig. 2a. The moment distribution in the elastic por- 
= tion of the plate can | be treated as a case of bending ofcircular plate by shear- 

ing force Q,. with superimposed bending moment M,;, uniformly distributed 


o along the inner boundary r=c as shown in Fig. 2c. _ The intensity of Qe is 
is found from the equilibrium condition and the intensity of Mrc is found from 
>. _the plastic solution for the portion of the plate within the circle | of radius “ae 
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— Elastic Regione 
— 

@ 

1 


For the plastic be c, , the moment can be found by 
the equation of equilibrium and the that is 


moment distribution is as 
(8a) 
At the plastic-elastic interface, r = ¢ and, aan = 


From ‘consideration of equilibrium, the shearing force Qe at r =c is found to 


27bQ) 


For the elastic region, r< a, the moment distribution can be 
7 mined by combining Eqs. 2 and 3 and is found that y 


an shearing force Qp increases. The moment distribution in the plate in both 
ac elastic and plastic regions is defined by Eqs. 8 and 11, and is plotted in Fig. 3 


a for " increasing ‘intensity of Qb. It can be seen that as the shearing force in- 
creases, the plastic expands from the inner boundary outward and there 
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CIRCULAR PLATES a 


is no point in the elastic region which the full plastic Thus, 

- we can conclude that, as the shearing force Qp increases in intensity, the 
plastic region expands from the inner boundary outward. When the shearing 

force reaches the limitin 
and the circular plate becomes fully plastic. 

Now let us examine the possible mode of deformation: Since the plastic 

flow may take place only after the material points reach the full plastic : 

dition,» there will, ‘be no deformation in the plastic region. . Since the stress 
= in the plastic region is in gee BC = shown in Fig. 1, the associated 


=0:1 


However, because « of the continuity seleiaas at the interface r = =c, no de- 
ae can take ‘place aa the plastic region as long as the plate is only ‘par-— 


The circular plate becomes completely when the force in- 


creases and reaches the 


as shown in Eq. 8 8. Then plastic plate is governed by E Eq. 


the distribution of deformation rateisfoundtobe 


Af(a-r)- 


Namely, the circular when the shearing force Qp reaches the 


2 


— — 
— 
_EM1__ 
— 
1 
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— 
j 
the simply supported circular | pate 


to uniformly distributed shearing seven Qp along the inner 


- According to elastic analysis, the permissible shearing force Qp which will - 
_ not produce a plastic region cannot be greater than: ae 


b/a = 0. 1/3, Eq. 16 becomes 


‘The carrying circular plates with various prescribed 
q boundary conditions are determined in the following by procedures described — 
in previous sections. Design charts are also constructed for various cases» 
¥ from the determined load carrying capacities. In these design | charts, the re- — 
- lations between the load parameter, p/p, the thickness-diameter ratio, h/a, 
and the geometrical parameter, b/a,are arranged so that the permissible load | 
can be found either from a given geometrical configuration or from the thick- . 


ness of the plate. is — in I to illustrate the use of 


Uni the Inner Edge Built in at 
Outer Edge.—Fig. 4 shows" a circular plate built- -in at the outer edge is sub- _ 
5 ba to shearing force Qp uniformly distributed along the inner edge. During © 
j 


plastic flow, the whole plate is in the regime C-D of the yield condition shown @ 


in Fig. 1. The stress state at the inner edge is represented by point C, and b 


the stress state at the outer edge is represented by | point D where the forma- _— 
tion of a plastic hinge is admissible. 4 Because of this, the plastic flow be- 
comes possible | when the plastic region reaches the outer r edge as the limiting ‘. 

load is plastic limit load is found to be: 


q 

7 % comparing Eqs. 15a and 17, it can be seen that the load carrying capacity a: Be 
[a  . of the aforementioned circular plates determined on the basis of limit design __ Sane 

[| qj is about three times as great as that determined on the basis of elastic analy- ae 

_CIRCULAR PLATES SUBJECTED TO AXIALLY SYMMETRICAL LOADING i! 
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of the plate is plotted in Fig. 
2 pate Loaded; Simply Supported at Outer Edge. —As oer in 1 Fig. 5, 


tie limit load is reached, the whole plate is in the ‘regime 1 B- C of the ad 4 
condition shown in Fig. 1. 3 The plastic limit load was f = determined by _ 


‘Eq. 19 is represented in Fig. 
Uniformly Loaded; Built- “in at the Outer Lies. shown in Fig. 
lar p'ate built-in at the outer edge is subjected to a uniformly distributed a 
of intensity p over a circular area of radius b. This case has been = enal 
gated by Hopkins and Prager. 4 When the limit load is reached, a portion “= “ils 


from r is tor= is in the regime B- Canda portion of the from 


in Fic. en For the case in which b/a< 0. 0.6066, the ratio c/o< 1.0 and the plas- os 


“tie limit load is" defined by the following relation: 


the case in in b/a< < 0. 6065, the ratio 0 1.0 the plastic 


ong pb 


; 


tion of the plate determined from 20 and 21 is in ‘Fig. 


4 _ When the plastic limit load is reached, the whole plate is in the ‘regime B- 4 of 
the condition shown in Fig. 1. The} plastic limit load is found to be: 
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is reached, the whol plate is in the regime -D of the yield condition 


 \a 


Eq. 23. is in Fig. 


il circular plate built-in at the inner edge and subjected toa shearing force ae 
uniformly distributed along» the outer edge. This case was first solved by 
Haythornthwaite. 7 Ww hen a critical value of the load is applied, the tangential a 
_ moment at the inner edge first reaches the magnitude of limit moment ae: 
3 However, it is possible for the plastic region to propagate from the inner edge 4 


4 nitude of the limit moment My as the load increases. When the limit load is . 
reached, the whole plate is in the region B-C of the _— condition a in 
1. | The Mastic limit load is found to 


relationship. given by Eq. 24 is plotted Fig. 9. 

ghee Force Along the Outer Edge; Built-in at the Inner ‘Eds re; Outer 

Edge Prevented from Rotating —Fig. 10 shows a circular plate built-in at the at 
_ inner edge and subjected to a shearing force Qa uniformly distributed along 
the outer edge which is prevented from rotating. | When the limit load is 
reached, a portion ¢ of f the plate from r jp to r=c is in the regime B-C and ea 
portion of the plate from r=c to as a is in the ; regime C=D of the yield con- : 
in Fig. 1 load is defined the following rela- _ 


_ Eq. 22 is represented in Fig.7 
— _ Circular Plate with Cut-out; Uniformly Loaded; Built-in at the OuterEdge.- (i 
Fig. 8 shows a circular plate with a circular cut-out built-in at the outer edge 
; 
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The r between the load parameter and the. geometrical con- 
fi figuration as defined by Eq. 25 is plotted in Fig. 10. eee ae eee i 
_ _Uniformly Loaded; Built-in at the Inner Edge; Simply Supported Along the 


edge is thn te from rotating but is free to deform. _ When the limit load is 7 
reached, the whole plate is in the regime B- Ce The plastic limit load is found 


to be: 


jected to to : a uniformly distributed load of intensity p. The outer edge is pre- 
¢ "vented from rotation. When the limit load is reached, a portion of the plate 
from * b to r=c is in the regime B-C anda portion of the plate from r r=c to x 
r“@ is in the regime C-D of the yield condition shown in Fig. 1. The radial 7 
moment is equal to zero at r=c which is defined - the nie ‘pen ae 
20 h 


e plastic limit load is related to the geometrical configuration a . 


— 
We in(2) = ln | 1 - 4 4 
| Outer Edge.—Fig. 11 shows a circular plate simply supported along the 
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it 
— 
— 
— «|. 


3 


OS 06 


< 
~ 


oO 


Janua 


1960 


af) 


Ooh 


or 
“Ae 


* 


— ae. ge 


By combining 27 and limit load can be and is plotted 
In the analysis described, the limit load of is deter- 
. AS ‘mined by taking only the bending moment into consideration. The material can | 
ss deform plastically and is assumed tobe rigid-perfectly plastic. Theoretically, — 
the circular plate collapses” when the plastic limit load is reached. For the - 
_ case in which the strength of the plate is of primary concern, the plastic limit a 
load so determined is generally higher and closer to the failure load than the 
allowable load determined by elastic analysis. Therefore, the use of the plas- 
ic limit load as a basis for the design of circular plates is more realistic and 
will improve the utilization of materials as well as the strength-weight-ratio, — 
= such an analysis gives the designer abetter estimate of the load carrying 
capacity of circular plates than the elastic analysis. 
Now the question is whether the plastic limit load of circular plate so de- 
termined is on the safe side. The answer is yes. It has been proved analyti-— 
_ Cally and experimentally that circular plates in reality can withstand loads — 
higher than the theoretical plastic limit determined in the foregoing. 
As pointed out by Onat and Haythornthwaite,8 > however, circular plates are 
2 strengthened by the monnrans stresses after they are deformed from their 
Initial flat configuration. Experimental evidence has been cited to support 
_ Furthermore, experimental results of R. M. Cooper and G. A. Shifrinl3 _ 
"show that the deflection of circular plates under the plastic limit load is of the 
same order of magnitude as elastic deformation. — This insures t that euenetre 
deformation is not a factor which need hamper the use of the plastic limit load © 
as a basis for — in cases where the strength of the cir cular plate is of 


paper presents design charts for subjected to axially 
ie loading based on the | plastic limit loads determined by the analy-— 
_ tical methed given by yng and Prager. 4 It was noted that the plastic Mimit~ 


Ras * Cooper, R. M. and G.A. Shifrin, “An Experiment on Circular Plates in the Plastic 7 
” Proc. 2nd U. S. National taal of — Mechanics, PP. 527-534, 1955. or 
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3 
Circular plate, in cases where the strength is of primary concern in design. 


5 a The use oft the design charts is very simple. One of the quantities a, b, 2! 2 by 
p, and To can be ‘determined after the others are given. After the boundary 
he conditions and the type of loading are known, a suitable chart can be chosen. 
For example, anozzle is connected to a pressure vessel as shown in Fig. 13. 
determination of maximum allowable internal pressure for this vessel 


000 x = = 1500 psi 


The he vessel can be asa circular in the inner 
“at 
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FIG. 13. 
‘tome to Fig. 13, the di- 
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© in the vessel is greater than that in the nozzle 
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edge subjected to a uniformly — 
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0.25. 


this of p/do, the plastic limit load for the end plate is: 


= 30,000 x 0.018 = 540 
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FOREWORD 
na This paper is the fifth in a series of reports | on bi cus 
Plastic Design emanating from a Joint Committee of ee 
Welding Research Council and the American Society 
of Civil Engineers. Reference should be made to 
Progress” Report No. 1 (which appeared in the July, 
1959, issue of the EM Journal of the ASCE) for such ee) 
‘items as the Guages, 
Nomenclature. e. For information, the complete ‘Table 
of Contents appears at the ‘beginning of this Progress 
Report. Chapter 7, entitled “Compression Members,” > 
‘The Joint Committee has approved the publication 
of this information with the strong recommendations 


Note. .— Discussion open until June 1, 1960. Separate discussions should be 
for the individual papers in this sy mposium. To extend the closing date one month, a 
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“a CHAPTER 7 


"Simple. pl plastic theory assumes that to 
ments will sustain a certain limiting bending _moment value (the “plastic 
“moment” Mp) that is dependent only on the geometrical properties of the = 
cross section and the yield stress of the steel. When this maximum moment — 
is approached, curvature increases indefinitely anda hinge type of action ocr 
curs, _ The presence of axial force tends to alter this situation. In the follow- | 7 
ing the various effects of axial force on the behavior of members ina rigid 
ee As loads are applied to a rigid frame, the individual members which com- 
prise the structure are subjected to various | combinations of axial thrust, end _ 
- moments, and end restraints. . In many cases in plastically designed struc-_ 
tures, the stiffnesses of the end restraints do not influence the of the 
column because of the formation of plastic hinges. mola 
For. an individual member which is subjected to bending moments and an 
axial force, and which is sufficiently braced in the lateral direction, or a 
_ member which is bent about its weak axis, the mode of failure as be that of A 
excessive bending in the plane of the applied moments 
‘ane _ If a member bent about the strong axis has insufficient lateral bracing, and 
ifa large difference exists between the bending stiffnesses about each of the | 
v4 _ principal axes of the cross section, the ‘member may bend out of the plane of 7 


the applied moments and twist at the same time. This type of failure is called 


de- 


id An additional type of failure, which is basically different from those 
she above, involves t the structure as a whole, and it may occur when cae 
‘Sway of the structure is not prevented. This condition is characterized by a 
«shift of the total structural deformation pattern from one that is symmetrical © 
to one that is anti- -symmetrical and is accompanied by an -over- -all iateral 
of the frame. ‘This ‘situation will occur in a symmetrical 
fe ture that is symmetrically loaded when, at a critical magnitude of the yg ae 
-— total resistance of the structure to lateral movement becomes zero. The 
possibility of this “side- -sway” or “frame” type of instability places a re- 


striction on the ranges of applicability of the individual member-strength 4 


further situation develop. when the structure in « question 
horizontally from the first load application. For these cases the horizontal 
displacement of the column top with ——- to the base may alter the carry- — 
ing capacity of the individual member. This type of action would be similar _ 

in nature to that limiting the maximum carrying capacity ‘of beam-columns. 


AS: in that case, ‘there would be reached a certain loading for which the struc-_ 


a ture would continue to deform in the direction of initial movement -(“exces- oe 


AS pointed out in previous sections of this Commentary (Art. 6.3), 2 
problem of rotation capacity may also influence the design. It may be mace 
sary in contain situations to insure a relatively large rotation at 

Tn the following articles, the effects” of axial force will be discussed, . 


will” their consideration in the design of one- and 
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member is to the combined action of moment and 
axial force, the” available plastic moment capacity is reduced from the full a 


value of M, to a lesser value that will be. designated as | Mpc. However, the 

design procedure may be easily modified to take this reduction into account 7 
“plastic hinge” charac- 

teristic is still retained at “this lower Because this is a property of 

_ the cross section, the value of Mpc as considered in this article is indepen-— 

dent of loading and the slenderness ratio, and it is immaterial q 
whether the axial force is in compression or tension. Even though instability 3 
ae are excluded, the methods developed here give a good approximation ; 
of the actual ‘behavior of very short compression members and 


other columns of practical proportions, as will be shown in Art. 7.3. 


Mie The first published work on the influence of axial thrust on the moment — 
capacity of a short column is that of Girkmann. (7. 1) Methods of determining _ 

_ Mpc are given by Baker, (7.2) Roderick,(7-3) Hendry, (7.4), Neal, (1. 6) and by pia 


Beedle, ‘Thurlimann, Ketter and Johnston. ( (3.1, 7.5) 


ots As an illustration of the influence of axial thrust on the plastic moment 
value, consider a rectangular section as shown in Fig. 7.1. Assuming, for 
example, that the thrust is maintained constant and that the moment is pro- 
4 gressively increased, the moment-curvature relationship will be that shown — 
; a non- dimensionally in Fig. 7.1. (The moment is non- -dimensionalized by di- - 
ne _ viding it by Mp, and the curvature is non- -dimensionalized by dividing it by By, aa 
; ‘the curvature at the inception of yielding. ) It is assumed here that the axial 
force is in compression. . For purposes of comparison, the curve for no axial 
=. is shown as a dashed line in Ae starts first in the outer 


As was the case when axial force was not present, the hinge eieniiahion cor- of 
responds to that situation in which the cross section is fully yielded. “(Stress = 
‘diagram D in Fig. 7.1, or as redrawn in Fig. 7.2). However, there will” 
7 longer be equal areas yielded in tension and compression as in the case of 
_ bending without axial load. Thus the neutral axis no longer coincides with the 3] 
_— centroidal axis. Hinge rotations for the computation of the virtual work per-— 
formed by the mechanism may still be assumed about the ¢ centroid of the sec- 
tion; the error for the structures considered inthis Commentary | is negligible . 
_ For very deep built-up sections, however, the effect of the shift of the neutral - 


of equilibrium of the ‘internal an¢ and external forces; that is, 


| | 
— 
— 
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In Eq. 7. 1 o is the stress at a given fiber, y is the distance of that fiber from 
centroidal axis, dA is the differential and the integration 


(see Fig. 1. 


* Yo! 
Substituting .2 into Eq. 7, 3, and noting that o =O0y and that 4 


no axial force present Mpc = Mp = = Oy b a?/4, | the non- -dimensional 


on 
) 
y’ 


$0 
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ait the moment increses: 
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By a similar process” the values of f Mp car can be computed for any 
; "Equations follow for wide- flange “sections subjected to bending about the | 
: strong and weak axes. The nomenclature used is that given in Progress Re. ” 
is port No. 1 “of the Commentary, and as shown on the inset of Be 1. 3 . New ; 
terms are listed in the nomenclature appended to this report. 


“Neutral axis axis in flange: 

A (1-P/P 


and Zz. 
_ Although Eqs. 7. 5 5 to to 7. 8 are given . in a forn form ‘suitable f¢ for determining Mp Po 
r= a given wide- -flange section, they. do not readily indicate the influence of — 
cross-sectional shape. _ Equations for strong axis bending = reflect this 
influence are given in Ref. 3. 1 in terms of of the ratios Ag/Aw- _ These equations ¥ 
- contain the simplifications that for wide- flange _ shapes d/dy a and ‘dt/dy are 


about the same for all used as and é are equal to 
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PLASTIC DESIGN 


Axial Load 
FIG. 7.2.—STRESS-PATTERN AT 


‘FIG, 7.3.—IN ~INTERACTION CURVE FOR STRONG AXIS BENDING 
WF SECTION (MOMENTS ABOUT X-X L= 
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in Fig. 4 The corresponding equations: for weak axis ending of 


Neutral axis in web: 


= 1. (w/b) (a 


0< 


(d/d =. P/P 1.00 


h typical values w/b = 04 


7“ 4 are obtained for weak axis bending. ae 


= 


9 


FIG. 7.4. INTERACTION | CURVE FOR WEAK AXIS B BENDING 
OF WF SECTION (MOMENTS ABOUT Y-Y ONLY; 
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the of formulating rules (see section on design : 
The behavior of wide- -flange columns subjected to. biaxial bending (that is, 
where the is not about the principal axes) is 


Fig. 7.5 shows the correlation between set of experimental data 
ee (5.2) Each of the three tests was carried out on 12 WF 36 ‘members _ 
which were so short that instability was no ) problem (L/ry = 7 1. 0). . Further- 


a0 more, end about the weak axis was prevented by the the ‘arrangement of the 


F IG. 7.5. _TESTS 
- imife edges at the column ends, This permitted rotation in the strong direc- 
“tion only. The results are plotted as moment versus curvature to indicate a 3 : 
- influence of axial force on the reduction of the plastic hinge moment. For 


- member T-1 the axial load was zero (pure bending). Members T-7 and T-8 © 


loaded by an eccentric force. As would be expected, the hinge moment 
T-7 and T-8 not reach the oe value of ‘The binge 


termined value of Mpc was close to that predicted by ‘Eqs. 1. 5 and 7.6 
sented by horizontal dashed lines in Fig. 7.5). 
: : To afford a clearer indication of the correlation between this set of test Se 
e data and predictions, Fig. 7.6 gives the experimentally determined initial — 
--yield and ultimate strength values plotted | on an interaction diagram similar 
to that described earlier.* Included in this comparison is the result of a pure f 


_ ¥The interaction diagram for the theoretical ultimate strength was constructed using 


Ss 5 and 7. 6; the initial yield curve represents the elastic limit computed from the ’ 


Ae 
— 
— p | | e510" 5-055 
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196 
axial load test (T- 6) of | the same cross "section. ‘Due to the presence of 


» the elastic interaction curve overestimates initial 


capacity; this is because of strain- hardening. 1 
Results of a series of tests conducted by Hendry(7. 4) on short 3 in. stand- 


ard British I-Beams are shown in Fig. 1.1 ; excellent agreement with the 7 
theory is indicated. A constant moment was first applied to the member “a : 

i axial force was subsequently increased until the member had fully plastified. _ 

ee _ The experimental evidence shows that for mild steel the reduction in mo- 
al capacity due to the presence of axial force can be closely salen by 


‘To account for the influence of axial force in ce any of the appropri- 
ate equations or. curves of this” section could be used. However, since 
Li curves for wide- flange shapes fall within a relatively narrow band (see Figs. 
4. 3 and 7.4), it is possible tor obtain simple approximate expressions for these ~ 
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A. If the axial force is larger, he 


“the corresponding for weak -axis bending is: 
tpe/Mp= 1.19 [1 - (P/P, A: 
12 need only. be used when P is more ‘than 40% of Py. These approxi - 
mations yield an error of less than 5% (Figs . 7.8 and 7.9). oa yk ster: UC}  &§ 
these interaction equations is to find a trial section 


filled. 


| 


Theoretical ‘Ultimate Load Curve 


fei 


COLUMNS LOADED BENDING A\ AND 


COMPRESSION. 
“Thus” a guide may be formulated in order | to account the 
ence of axial thrust. It must be kept in mind that instability effects are 
ae neglected, and so these equations apply in the strictest sense only to short a 
columns. Calculations including the effects of instability have shown that the 
~? strong axis interaction equation (Eq. 7.11) is also valid for columns subjected =, 
‘In summary, the following equations are suitable for determining the re- = 
= duction of moment to axial force: 


— i 
rrections until the conditions of the interaction equations are 
| 


REDUCTION OF PLASTIC MOMENT CAPACITY 


For 0 <P <0.15 P,, 


For 0.15 Py <P <P 

A Weak Axis Bending, WF 


For 0 <P <04P, » Mpc = 
| For 04Py <P < Py, Mye = =1.19 | 1- - (P/Py)? | M, 
Mpc = | 1-(P/P 


re 
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1. 3 MOMENT- CARRYING CAPACITY 0 OF COLUMNS 


Although the solutions that were in n the preceding section repre- 
os a basic characteristic of the cross section and are adequate for pn 

columns, they do not always correspond to the amount of loading that a — y 
As an illustration of the problem, consider the eccentrically loaded cohima 
of Fig. 7.10. As the load is increased beyond initial yielding at et or 
_ plastification progresses along and across the column, thereby reducing its 7 
resistance to further loading. Fig. 7.11 shows a typical load-versus-center- _ 
ea a curve for an eccentrically loaded column. The portion of the curve — 


th 


WIDE- FLANGE SECTION (WEAK AXIS BENDING ON FC 
from OtoA represents the column behavior when the stresses are still 
he The portion A - B represents the range of partial yielding. — Finally, when 


| 


the P-versus- -yg, curve reaches point B, a further increase in load becomes 
impossible because the internal stiffness of the column is just enough to re- 

— _ sist P and the moment P (e+yg, > it is with the determination of one's maxi - 
‘mum value of P that this section of the Commentary is | concerned. 

% dent that this type of failure occurs by virtue of excessive bendingin the plane 


ue column will not fail i in the manner just discussed i if the cross section. 


&§ 
— 

— 
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two pirncipal directions (a characteristic of wide - sections). If the 
member is ‘subjected to bending about the stronger of the two axes, _and if no ~- 


_ strength will be reduced. This phenomenon is known as | “lateral - torsional 


uckling” and is the subject of Art. 7.5o0fthis chapter, 
2 If lateral-torsional buckling is prevented, the eccentrically loaded column Ls 
or the column with end bending moments plus an axial thrust) can “fail” only 
after a ‘certain amount of yielding has taken place. The problem is therefore 
not one of stress but one of stability. point of indifferent equilibrium | 
curs when the internal stiffness is just enough to resist the external moments. : 
= for certain columns the condition of = plasticity (defined by stress te 
am D in Fig. 7.1) will notbe reached. 
ce) solve the problem of the eccentrically 1 ery a 
on the basis clanateal stability was von Karman. Bleich has pre- 
ie sented an extensive résumé of the early work on this topic, , including the work oe 
of von Karman, Chwalla, Westergaard, Osgood and Jezek; it is contained in 
sections 9 12 of Ref. 6. 


column will twist and bend out of the plane of. loading; in 


LOADED COLUMN = ECCENTRICA LLY LOADED COLUMN. 
ot Recent work has extended the investigations of the early researchers i 7 
_ this field. Solutions to the problem of eccentrically loaded, end- restrained — : 
rectangular columns were presented by Baker, Horne and Heyman(1. 1) and by 
Bijlaard, Fisher and Winter. (7.8) Ketter, Kaminsky, Beedle and Galambos 
investigated the problem of hinged-end columns of as-rolled wide -flange 


| bent about their major axes.(7.9, 7.10, 7-11) 


ioe Parallel to the efforts of the determination of the strength of Be > 


‘loaded columns by rational | "means, attempts to define empirical 
_ equations for predicting column strength have been made. Among these , the : 
most recent and most extensive work that of Campus and (6.32) 


equ 


a 
a. 
ff 
, 
= 
— 
— 
| ~ Mtl ~ P/ < ‘1.000 


this value includes the influence of the effective length. 
Mequ = equivalent end bending moment which may be expressed as follows: 


larger of the two end moments. 


“Te smaller of the two end moments. T The sign of M2 is negative in 


Eq. 7. 14 if it causes a curvature opposite to that caused by M)- 
= the elastic > (“Euler” ) buckling load of the column in the plane of the — P 


moments 


full plastic moment of the cross section. 


Ye In case the maximum ‘moment occurs at the end of the. ec the f follow- 2 % 


"notably by B Bijlaard, Mason, Fisher, Winter, 8, 7. 


Johnston, Cheney, Ketter, Beedle, Kaminsky, (7.5, 7.9, 7.13) Campus, 
sonnet (6. 32) and others. (7.14) It is impossible in an article such as this to do 


5s full justice to the immense effort that has been put forth in n solving the v te a 


aspects of the eccentrically loaded column problem. Therefore, only a few | i 

have been mentioned; for a more complete review of all work, Ref. 6.19 lists en 
work up to about +1950 and a forthcoming publication Column 
search Council, (7.15) will summarize the most recent research. 


From the many available solutions, one has been selected which ‘stguesnan - 

‘most nearly the conditions that exist ina plastically designed rigid frame of f 

the type considered in this Commentary (that is, one-or two-story frames). 

_ The results of Ref. 7.10 have been chosen, because they represent a so- -called © 
“exact” solution (for including the influence of residual stress on = 

strength) and are directly applicable to columns fabricated from as- rolled 


a a) The mode of belies will be that of. excessive bending in the p af of t the 7 4 
moments; this plane is taken to be the strong | 


mi The material is mild, structural grade steel such as ASTM- et or A373, 


which is to possess the idealized stress-strain curve of Fig. 


an Members are originally ‘straight, free from accidental end eccentrici- ; 


Mec 
& 
— 
3 
| 
ak 
— 
— 
— 
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e) before bending remain plane a after 
f) The end slope of the deflected column is small. — 


a: The vahnstiod’ of | the column in the frame is the same as that of an iso- 


- lated member loaded with axial force and end bending moments. | 


Interaction curves for strong-axis bending of a rolled wide-flange section 
are developed in Ref. 7.10 by an iterative procedure. The influence of re-_ 
 sidual stresses due to cooling after rolling is included in these calculations, .. 
solutions for as-delivered columns. particular “residual 


7.12) is typical of commonly used wide- 


y is assumed in n the calculations. 


_ Interaction curves for wide - -flange evmtene: are shown in Figs. 7. 
7.14 for two loading conditions. These curves have been computed numerically | = 
_ for the 8 WF 31 shape; however, they represent a good approximation for any 
a other wide- -flange section. “The | 8 WF 31 section has been chosen because of 
sits low shape factor (f 10 as compared to the | average value of i. 14). 
Hence: the curves are conservative other sections in the ratio of 
shape factor to 1.10. The interaction curves show a non-dimensional plot of a 
the between the axial (abscissa) and end bending 


stress of 0.3 0 


FIG. 12. -— ASSUMED COOLING ‘RESIDUAL 


(ordinate) f¢ 


onstant values of the slenderness ratio in the direction of 

re. TP 13 gives the curves for a loading condition in which two equal end 
Moments cause the column to bend in single curvature; Fig. 7.14 shows the 


complete description of the necessary to obtain these in- 
_ teraction curves is given in Ref. 7.10. Only a brief outline of the procedure =a 
is given here to show how one point on a curve is obtained. For a given sec- Ee 
YP. tion an axial force, length and end bending moment value are assumed. _ The a 
a end slope corresponding to this loading is computed by a numerical integra- ae 
Ee tion process; that is, curvatures obtained from the moment curvature diagram 
of Fig. 7.15 are integrated to . give the deflected shape and thus the end slope 
of the column. This process is repeated for several values of the end mo oo 4 
; ment (length and axial force remaining constant) until an end moment versus 
end rotation curve can be constructed. ‘The maximum point on this curve se 
corresponds to the highest end moment which this column can support, thus 
giving one point on the interaction curves of Figs. 7.130r 7.14.00 2 
Eioat The curves of Figs. 7.13 and 7.14 have been computed for steel with a yield © 
stress of 33 ksi. If itis desired t to use these curves for steel with a different a 
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— 
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ae stress level, an approximate result can be obtained if the actual ; 


slenderness ratio of the column is multiplied by the factor where oy is 
stress of the particular steel (expressed in ksi). 


_ Figs. 7.16 through 7.20 show the experimental correlation of various tests - 4 


‘with the interaction curves of Figs. 7.13 and 7.14. 


i Fig. 7.16 correlates the theory pith tests pe performed by Mason, Fisher and 
- Winter on hat- shaped sections. (7.12) This cross section conforms most near-_ 
ly to the assumptions made in the derivation of the interaction curves fae 
= ly, that lateral buckling cannot take place since bending is about the ‘week 


“ING OF WIDE-FLANGE SECTIONS 


welll asin the following figures, the “continuous line ‘refers’ to the ‘theoretical 

“curves for the given n eccentricities; the | points represented by circles, squares, o 

or triangles are the experimental points. These have been corrected to in- — 

clude the influence of specimen yield stress other than 33 ksi. In Fig. 7.16 — 

the experimental points fall slightly above the theoretical curves; thisisas 


n in Fig. 6 of Ref. 7.12. 
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would: be expected since factors of the hat sections are 
above those of wide- flange sections, and the residual stress is somewhat be- a 


na low that assumed in the theory outlined here (f = 1.18, 1.25, and 1.17 for the 


sections that were tested; Sue ( oy y= 0.23 oy by me measurement as ‘compared — 
0.3 oy 


__ The correlation with tests conducted by Johnston and Cheney(7.13) is shown 


in Fig. 7.17. The test setup was such that the columns were essentially pin 
er. with respect to bending in the strong ‘direction, and fixed ended in the 


/ 


4 


= direction. This v was done by | the use of knife — placed perpendicular — 
It 
ne is interesting to note however, that except for the tests which fall close to the = _ 
= where failure would | have been due to “Euler buckling” in the weak a 


direction (see dotted curve) the correlation with the theory which neglects 


lateral- torsional behavior is ‘reasonably good. constructing dotted 


> 


| 
Mective tenet W h was deter 
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predictions in ‘Fig . 7.18. a The DIE peeliien, As which the test columns were | 
made, are geometrically ‘similar to American wide- flange profiles. The end " 
conditions of the columns were essentially pin ended in both principal direc- 
tions, since the end fixtures consisted of almost frictionless, hydraulically- 
: seated steel hemispheres. For such end conditions the lowest possible re- 
straint is offered to lateral-torsional buckling. As shown in Fig. 7.18, most 

of the test points agree rather well with the theory which neglects this type ot 
buckling, even though failure was by lateral-torsional instability. —Compari- 


son was made - only for the seated case where 1 moment was applied at one end 


— 10 


Tat. —COMPARISON OF COLU MN WITH THEORY. 


for the tests: of Johnston on anh Correlation with: the 


indicated i1 in this” figure. lateral braces 1 were to 
buckling (except for tests —* “23, and T- 


ii 
Tests ion. The test results are compared w 
— 


of these braces was determined in accordance with Eq. 6.25 


Progress Report 4 of the Commentary. _ These braces were adequate for the 
prevention of lateral-torsional buckling. Good correlation is seen to exist 
_ between theory and experiment | in Fig. 7.20, except for test T- 13, which was was 
considerably stronger than predicted. "This was a relatively short column 
with a low axial force, and therefore its strength —— reflect the influence _ 
The above comparisons of theory with results show that the 
; = interaction curves of Figs. 7. 13 and 7.14 are capable of predicting column be-_ es 


» 
EULER BUCKLING 


WEAK 


FIG. 7. 18. —COMPARISON OF COLUMN TEST RESULTS WITH THEORY. Rav 14h 
cimens used the experiments. in all one of the test programs were 
rolled wide-flange shapes, probably containing a similar residual stress pat-_ 
. tern to that present i in the 8 WF 31 beng used in deriving these curves. The 


of Ref. 7. 12). 


ee 
| 
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design. As a convenience in interpolating, the interaction curves have been 

_ reduced to approximate design equations by means of curve fitting.(7.10) In 

order to arrive at simple approximate equations, their application is limited 

f tod< P/Py < 0.6 andto0 < L/r < 120. These limits place all practical — 
pes problems in plastic design well within the scope of the equations, at leas! for 

buildings oneortwo storiesinheight. 

In order to apply the following design rules, the column must be sto ; 
quately in the ‘Grectios to ‘the momenta so as to to 


Hn 


lateral buckling. The ends are assumed torsionally 1 re- 


j erally the case. For weak axis bending the use of the suggested approximate ard. 

formulas would lead to conservative designs. In this case lateral- iis 
¥ torsional buckling would not be a problem. | The above defined limits of axial _ 
ee force and slenderness ratio are appropriate, provided the whole frame is re- 


support of an adjoining structure of known stability, or ape stiffness 
furnished by the and the : roof. 
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FIG. 1. OF COLUMN TEST RESULTS WITH THEORY. 
pgoing solutions of the case where | the moment ratio is +1 (Fig. 
7.13) and 0 (Fig. 7.14) * may be extended to cover all cases of column loading 
inthe following manner: 
A column which is bent into double curvature (into an S- -shape) by 
a end moments (moment ratio equal to -1), can be thought of as two columns of 

_ the type shown on the inset of Fig. 7. 14, where the inflection point at the cen 7 


provides the moment-free pin Thus, if the slenderness ratio of the 


— 
— 
iia 

— 
_ 
 iargerone. It is positive when each moment causes the same curvature (single 
ae — wy. ture) and negative when one moment causes opposite curvature to the other (double ie ee 


‘column is | 80, the corresponding ‘slenderness of the effective column is only” 


m1 Oto 40 and for 0 <P <0.6 Py the curves }are almost identical to the curve 
= for full plasticity (L/r = 0). Thus, in the case of loading producing double a 
curvature, ‘the equation for a zero length member may be used (Eq. 7.11) 
without causing an appreciable error. . Under certain conditions however, this 
column may “unwind” from its antisymmetrical deflection configuration and 
i fail by concentric buckling. This occurs when the axial force equals the criti- J 
cal buckling | load for the corresponding axially loaded column. ‘The strength 
of the column ag against this type of failure is reduced due to yielding at the © 
ends where the moment is applied, and thus “unwinding”(7. 17) will occur be- 
fore the “Euler” buckling load is reached. To safeguard against this type of 
the application of Eq. 11 is limited to 0 P/Py <.0.6 
For columns with a moment ratio from - 1.0 to 0, the equations for the — 
more severe case could be used (that is, for a moment ratio ¢ of 0). Similarly, 
-— > equation for a moment ratio of +1.0 could be used if the ratio ranges from 
summary, the following guides are listed for t the design of steel 
a. columns | bent about the ‘major axis in frames where sidesway is pre- 


FLANGE COLUMNS, ‘STRONG A AXIS. BENDING 
= 


The slenderness ratio with to the axis of be 


to or less than 120 and the axial force must be equal to or less ide = 


0. 6 Py. | The column must be adequately supported laterally, and sid 


sway must be by. external 


Equations: 


1. 18 a- P/Py) Mo/Mip | 1.0 
Ratio of end moments from - -1 to (Fig. 7.14) 


Mo/ 


B- G 3 (P/Py), <1 1.0 


1.00 - K(P/Py) - -5(P/Py)2 
K = 0.4204(L/r) /70 (L/r)? /29, 000+(L/s)3 /1,160 1, 160, ,000 
7s 0. + 700- 00-(L/r)8 /606,0 00 


imi 


aj 


a 40. In Fig. 7.14 it may be observed that for slenderness ratios ranging a { 


, 
— 
‘a af i | 
44 
— 
— 
; 
— 
— 
(7.16) || 
— _G= 1.11 +(L/r) /190 -(L/r)2/9, 
— 
— 
— 
— of P and L/r should lie between 0< P<0.6Py — 


ae Eqs. . 7.16 and 7.17 ‘the quantity 1 Mo should be taken as the larger of the eS 

end moments. Values of B, are tabulated for slenderness~ 

ratios of 0 to 120 in Ref. 1.3. 5 


be ee. undergo a certain amount of inelastic ‘rotation; at the same Be 

time the axial load and end moments must be maintained on this column. This — 

— especially true of a column in which contains the first plastic hinge, since [i 

a this hinge must continue to rotate until the last plastic hinge of the frame is a 
developed. The problem is thus to determine whether a column 


a =z. needed rotation and me. ill support the loads for which it was de- 


rotation hinge in order to form a 

_ mechanism under a given loading condition can be determined for any frie e 

column by elastic-plastic method which has been developed by Driscoll. 18) 

- method will be discussed further in Chapter 9 of this Commentary. - 

a _ The theory outlined in Art. 7.3 of this chapter permits the calculation a 

the end moment versus end slope curve up to the point of maximum moment. | sa fe 

Information from a study currently underway(7- 19) ) will ‘permit the construc- 

_ tion of the full end moment versus end slope curves for any rolled wide- -flange . 

section. A theoretical evaluation of the rotation capacity problem v will there-_ 


evaluation of the performance c of columns in tests shows that adequate 
“rotation capacity will be exhibited for the usual one- or two-story structure, 


"provided that lateral-torsional buckling of the columns is prevented. __ 


facts may be deduced about rotation capacity from 
1) In all of the experiments on full-scale frames, described in Art. 5.3 of the ‘i 


_ Commentary, the maximum n loads which were sustained w were > at least = aa 


and dimensions of these frames were ‘of proportions. 
2) Observations of the rotational behavior of as- -rolled ‘columns 
: subjected to bending moments about the strong axis were made in Ref. 7.20 pe J 
ns _ The tests were performed on 8 WF 31 and 4 WF 13 columns, which ranged — 


in slenderness ratio from 27 to 111. Most of the columns were subjected — 


toa constant a load of 0.12 Py and to a number of different end mo- sf 
-ment conditions. No lateral was provided between the supports for 


of the “tests; the prevalent mode of failure was therefore lateral- 
torsional buckling. Many of the tested columns exhibited a rotation 


‘moment versus end rotation curve was still rising when the test was discon - me 

3) ‘The results from one of the current series of tests at Lehigh University 
are shown in Fig. 7.21. Also shown are the data pertaining to this test. 

end moment versus end rotation curve of this test is shown. 
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M 
. 6.3 of this 
6. The total rotation was about three times the rotation at 
- elastic limit when the test was discontinued. Even though the slenderness 
ratio and the axial force were considerably higher than those commonly a 
_ encountered in the types of structures considered in this Commentary, _ a 
end rotation capacity of this column appears satisfactory. Further tests 
on laterally restrained columns have shown that rotation capacity increases © 
d slenderness ratio decrease. 


getation For one- and two-story rigid frames the column ends can 

LS expected to supply the required rotation, provided lateral- torsional | buck- = 
ling is prevented (see recommendations in Art. 7.5). For the very infrequent ee 
cases when the axial force is high (P > 0.3 Py), when the column is very = 


slender (L/r A 100), © or when the equal ¢ or nearly equal end moments cause 
single cu 


if it were found that this was not the case, one could secity a larger c column — 
in order to force the to occur elsewhere first. 


1.5 THE INFLUENCE OF LATERAL- TORSIONAL BUCKLING 


Ifa column has sufficiently different bending stiffnesses: in its two princi- 
_ pal directions and if the external bending moments are . applied in the stronger re j 
_ direction, the column may not reach the strength implied in Art. 7.3 unless 
adequate lateral bracing is provided. It will usually fail by lateral- -torsional — 
‘buckling before excessive bending in the plane of the moments is reached. ates 
‘This type of buckling has been observed by investigators who have con 


5 


‘one - or two-story rigid frames, lateral buckling does not take 


ing. a a inelastic solution has been obtained for a column bent by equal end 

moments which cause single curvature deformation. (7. 23) Good correlation - : 
of this theory with test results was observed. Further ‘research is 

‘directed towards the solution of other loading cases. 

_ Lateral-torsional buckling of a column which is located ina 


The rotation n capacity may be 


*For a summary of these see Chapters 3 and 4 of Ref. 6.19 and Ref. 


— 
& 
— 
— 
| 
— 
— 
4 
— 
| 
im 
elastic lateral-torsional buckling are available for a wide range of loading 
— 
J 
= structure may tend to exaggerate the following two effects; a 
ory of Art. 7.3 may not be 
4 1e greatest reduction in strength is associated with columns deflectedin _ 
single curvature by equal end moments if the ends of the columns are — 
_ pletely unrestrained in the weak direction. (7.23) Tests of columns under 
ee! @é€8§©§=6equal end moments causing single curvature but in which almost full re- ra al 
peo 


the axis was achieved, (7. 13) ‘show that in all, in- | 
stances it is possible .to reach the strength predicted by the theory of Art. 7. 3 7 
: (see Figs. 7.17 and 7.19). The exception is where the axial loadis close to =. % 
- the “Euler” buckling load. For other loading cases such as the one shownin | 


‘Fig. 7.18, , the reduction of strength due to lateral-torsional buckling i 


_ It may thus be concluded that for the type of columns which occur in the 
structures considered in this Commentary (that is, columns for which the 
axial force is relatively low), the column theory of Art. 7. 3 will give a satis- > 

of column strength provided that considerable restraint of 
the column ends in the weak direction is provided. | - Commonly used base- — 
plate and anchor -bolt arrangements, and the necessary transverse beams be- 

— adjacent frames at the column top will in general insure this. sae) Mies 

_ The second possible effect of lateral-torsional buckling is to influence the 
setntien capacity. Since Ref. 7.20 shows that the rotation capacity may be re- 
duced if the ‘member fails by la lateral instability, it ‘recommended that 
lateral-torsional buckling be prevented by suitable intermediate bracing in 
regions of plastic bending. _ Until future research discloses a perhaps more — 
liberal procedure, it is suggested that the spacing ef the bracing points be = 4 : 

- eaeaie with the recommendations made in Art. 6.3 (Eq. 6. 25). In case 
= loading produces symmetrical single curvature, the minimum bracing ce 
tance of 35ry s should be used throughout the column length. Fig. 7.21 shows — 

the results of atest in which the lateral support spacing was in accordance © a 
with Eq. 6. 25. — was to lateral “torsional 


i cobalt subjected to an axial force and to bending moments in one of e 
principal planes can be predicted if no lateral- torsional buckling occurs. 


can carry may be less than the ae computed on the basis of the be tr = — 
_ individual members if sidesway | is not prevented. In this case the a 
bility” exists that the frame as a whole becomes unstable before the plastic 

is formed. If this occurs the structure te said to have failed 


a If a frame is combination of vertical 
forces, it deforms laterally from the first load application. The change ec 
in geometry may alter the carrying capacity of the individual column, a 
since the column top is no longer over the column base and hence ad- es 


ditional bending moments are introduced by the vertical loads. 2 a 


whole » structure becomes unstable in this deformed | position much like 
the loaded beam-column of Art. 7.3. A load-deformation 


‘ 


28 
-*If a calculation of the reduction in strength is desired, the methods of Ref. 7.23 may ¥ : 
be ‘used. . An 4 ee of this strength may also be obtained by using the interaction _ 
p is nen by the oriios’. lateral buckling moment, and 
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1960 
curve i is ; shown for this type of instability 1 in Fig. 


critical loading, the structure continues to deform without a an increase 


wre 


9) If a symmetrical st structure is loaded by , symmetrical vertical forces _ a 


hs. 7.22), it is possible that the frame may pass from a symmetrical, Fi a 
stable deflection configuration to an unsymmetrical, unstable con- 
om figuration. At this instant the total resistance to an over-all lateral 
movement becomes zero. The load- “deflection on curve is characterized 


he a shift from a a situation where P can increase as the deformations an 


crease to one where large deflections develop without an increase al 


This behavior is analogous to that of a centrally loaded column, 


ae in which bifurcation of equilibrium is possible at a critical iar 


Ma faximum n Momen =0.661M 
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FIG. a. 21, MOMENT VERSUS END-ROTATION OF OF TEST A- 


cae “The first of these | phenomena is now under study by Merchant.(7.24) These 

- theoretical and experimental studies are attempting to establish | limits of _ 

a loading and frame geometry in which frame instability will not be a problem. i 

_ ‘They will have particular application to multi-story frames. This type of 
instability has never been observed in any of the previous full-scale frame 


floor slabs and the roof. None of the aforementioned | full scale frame — — 
received benefit of this additional stiffening. Nevertheless, they pe 


| 


PLASTIC D DESIGN le, 
second effects under extensive investigation 
~ structures where the forces may be assumed to act in the column see 


when primary bending effects: can be neglected (Fig. 7. A summary 


4. 25. In such cases the “effective” length of the column is said to gol 
creased beyond the actual length. it is possible for the effec-_ 
__ tive length to vary from |} kL= L/2 to k | L= . Usually the restraint offered — " 
the effect of walls, roofs, and floor is sufficient to prevent frame 
-—metabaltys the forces necessary to prevent it being relatively small.(7. 26) 
s The most serious instability condition results when the column bases are 
z _ actually pin ended. ‘It should be pointed out that partial base fixity, a oe 


Point of instability 


re 


axial loads in the columns and little or no bending) is seldom met in a 
- Rigid frames are primarily designed to support loads by bending action ss 
rather than by compression. Chwalla has shown (as discussed on p. 225, Ref. ‘i 
6.19) that for a symmetric, single-story, single-bay rigid frame, which ao is 

| es in the beam by two equal, vertical forces at equal distances from a 


load obtained by | placing all the loads on n the columns, despite the presence of ; 
i. _ At the present time there is no method by which inelastic frame instability — 
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‘Plastic Real Hinge- 


elle 


>. 
ne 


Buckling Ae 


f 


methods will soon. "Experiments. are being 


dicted on steel rigid frame models at Lehigh University to determine the . 
iste "range of the axial force and the slenderness ratio for which frame instability 


Approximate calculations show that frame instability for plastically de- 
4 signed structures will not occur if the slenderness ratio and the axial load > a 
ratio of the columns are held within the following limits: 
Justification for ‘this rule ‘may be found by -story, § 


-bay shown in Fig. ' 7. 23. Just before takes piace, par- 


— — 
— 


tially plastified zones ms at both knees of of the fi frame. At the 
’ “windward” knee unloads (that is, it becomes elastic again) and the “leeward” 
knee becomes more plastified due | to a corresponding increase in rotation. 
; ‘The curves in- Fig. 7.23 show a safe solution of this problem wherein this” 

plastic hinge is replaced by a “real” ” hinge (resisting moment neglected) and — 
the loads are moved to the top of the columns. . The straight line of Eq. 7.19 


is seen to be safe when compared to the curves for the analogous frame for _ 


the ratio of Ic h/Ip L < 2.5. Partial base fixity would increase the — 
ance of the frame considerably as shown recently by experiments on frame | 
models at Lehigh University, where partial base fixity was by 
rigidly attaching a beam between the ends of the column 
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= Effective column length coefficient 3 
Beam-column interaction equation 
1oment on beam-column 
ae Larger of the two end moments ona beam- POR adh 
= Smaller of the two end moments ona 


Critical load on axially loaded column 
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Closure P. Jr. and Herbert Ire 


4 ASCE.—Mr. “‘Augusti used infloence lines and integration ‘the 
“nt "error ¢ of the approximation in assuming that a distributed plastic curvature is - 
concentrated at the point of maximum moment. Actually, as stated inthe a 
for any given span, including non-prismatic and curved spans, the location o 
_ this concentration for exact results is at the centroid of the distributed curva _ 
ture. This principle is based on the linear relationship between the location n 
4 any angle change and the end- moment or force caused by the angle change. i 


This for the general case by the first terms of 


must always be linear. 
Therefore, the error discussed by Mr. " Augusti is simply the error in as- _ 
suming the centroid of curvature at the point of maximum moment, and it can — 


Usually be dismissed by inspection. As the paper stated, to eliminate 


wat tolerated for the illustrative ancegien. the authors would certainly agree with Ri 
Mr. Augusti that this error may occasionally be significant. 

_ Mr. Renzulli’s analytical determination of the M - - grelationship for a pa 
tangular section of an ideal material, being outside the scope of the Py will 
not be commented on here. _ References were listed in the paper for work in 
this area for steel and reinforced ‘Mr. Renzulli’s expressions for 


ing: asa mathematical expression of the problem to be solved. To be. cael ' 
ly general, terms expressing the joint angles from sway should be added to 

this equation. As the discusser stated, direct solution of the system of non- Fy 

linear simultaneous equations obtained by writing this equation for all joints a a 

practically impossible. _ He therefore recommends that the non-linear ‘terms 
= a estimated using an assumed set of moments, and that the resulting linear 4 
B . system of equations then be solved for the moments and sways. _ Generally, | the 


‘moments and s sways so obtained will not agree. with the assumed moments and s 
. sways, and new trials based on new sets of assumed moments will be neces- _ 
» sary. Actually, this algebraic procedure is equivalent to the numerical pro- as 


Dass of the paper, the only important difference being that in the latter the ¥ 


i 


December, 1958, by Lawrence Johnson, Jr. and Herbert A, Sawyer, Jr. a 
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moments of each trial are obtained by a ‘moment distribution rather than by a a 
olution of simultaneous equations. The superiority of moment distribution ces . 
an algebraic method is well known (although with electronic computers 
this reduced). For example, for the algebraic s solution of Ex- 


whereas the corresponding solution by moment distribution occupies four 
ae On the other hand, the interesting alternative approach to the problem sug- 
gested by Dr. Rosenblueth, involving assumption of elastic hinges at the plastic 

_ regions, is like the methodof the paper in that it uses moment distribution for 
= trial. To the authors, however, an important point of difference is that 


tong 


‘ae “over factors of all hinged spans for each trial. These elastic ‘quantities r 
‘main constant for the method of the paper. f rf fren 


_ The design of a structure usually involves analyses for several ‘tending 
conditions, and several moment- distributions are required for each elasti- 
analysis. Because of this repetition, the authors recommended and 
- illustrated use of unit-moment distribution factors. as the most efficient method - 
of moment distribution. The authors agree with Dr. Rosenblueth that the ad- yi 
vantages of this method diminish as the number of joints increase. _ However, *y 
this method is and the usual method of ‘moment distribution be 


_ deformations are neglected, Fig. 8(b), and the latter a method following Baker ae 
a in which thereis a sudden and complete transition from elasticity to plasticity — 
- without strainhardening, | Fig. 8(c). . Both of these methods are compromises = 
, the general elasti- -plastic method of the paper and the well-known | 
ra simple rigid-plastic theory, and the potential value of these compromises is — 
: that they may, unlike simple plastic theory, account for deformations, and that — 
- this accounting may be accurate enough for design purposes wi 
ee _ The authors cannot agree that Mr. Renzulli’s approximation involving neg- — 
; lect of elastic strains is generally good for “homogeneous and reinforced con- 


crete sections”. The ‘ratio of total elastic to total plastic. angle change for a < 


member with a linear moment diagram is equal to the ratio of the elastic to ie 


plastic areas below t the level of M of the M - -¢diagram, authors’ Fig. Ma) .. 
For the M - -> diagram derived and presented by Mr. Renzulli in his Eqs. (1) 
pt (2) and his Fig. 2, using his “ideal” stress-strain relationship with i : 
Sek this ratio has been calculated to ‘be about 9 for M= My, and ‘it ap- . 


As an example, for the M - -¢ relationship - 
ie = in Example 2, based on tests of reinforced concrete, the minimum _— 


‘ ‘the approximation of Fig. a I should not be used for reinforced concrete. ae 
As for steel, tests recently conducted by one of the authors at the University a 
of Connecticut yielded, for rolled steel sections with linear moment diagrams, 
‘minimum ratios of elastic. angle to plastic angle of from 0.09 to 0.70. _ (Since 
‘sections with rather slender elements of these the 


‘ 

3 
— 

a 

| 
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lower of these ratios is s probably more of the ‘usually- 
sections than the higher.) | Although ‘these ratios are much lower than those 
for reinforced concrete, other considerations cause the authors to wa | 
the use of the neglect-of- elastic-strain approximation for steel. First, an im- 7 
portant purpose of a strain analysis is presumably to detect those structures 
for which a limiting strain may prevent attainment of the ultimate loadindi- 
_ cated by the rigid-plastic theory. Unfortunately, those members with the small- a 
est limiting strain are also the members with the nighest ratios of elastic to 
plastic angles. Hence, the error of this approximation is greatest (and on ‘ti 
_ unsafe side) for the very members or structures for which a strain analysis — 
may be important. A second, and secondary consideration is that the | ‘post- — ¥e 
elastic portion of the moment curvature for rolled sections, 


not be accurately represented by any single line. 
_ The two other possible purposes of a strain analysis of a a steel structure in “* a Bi 
the plastic range should also be mentioned: one is the determination of the 
- amount that strain-hardening will allow the load on a structure to exceed the 7 
-fully- plastic ultimate load, and = second is the determination of deflections. 


ae 


Baker 


clined strain-hardening or “modified plastic” line, as shown in 1 Fig. 1(b), a =e 


FIG. 8. —COMPARISON OF M - CURVES 


For both of these purposes, the above errors _ would make the use ” the ot 
neglect-of-elastic-strain approximation for steel v very 4 
me The Baker-Rosenblueth approximation avoids the above errors by account- 
a : ing for elastic strains, and thus the authors would agree that this approximation — 

_ is of potential value for reinforced concrete design. . Since one of the authors _ a 
has recently discussed the relationship between this approximation and the a 
aa -elasti-plastic methods in these pages, (15) further discussion will be limited to — 
two observations. One is that the strain- -hardening portion of the M -@ ae | 

E lationship, neglected in this approximation, is essential for any rational de bi * 
termination of the plastic - region curvatures at the time the limiting strain is 
attained, and therefore essential for rational determination of ultimate hinge-_ i 
rotation capacities. (2) Second, the authors believe that this advantage and the \s 
other advantages(15) of elasti-plastic methods which result from use of the ful 
_M - ¢relationship do not seem to be outweighed by any increase in computation-_ 

al difficulty. It is primarily a lack of definitive information on M -¢ a: 
_ ships which presently ‘discourages use of elasti-plastic methodsfor reinforced _ 
‘The authors appreciate Dr. _ Rosenblueth’ ~ elaboration on the possibilities, 

_ Valuable in design, of at intermediate stages 
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‘January, 
= 
‘limit Mr. mentioned, the elastic properties used with 
ee plastic method are independent of the magnitude, location, or quantity © 
ao of the plastic regions, so the designer can make adjustments in Finforcing 


stool without and of procedure 


_ distribution, with its attendent advantage of pictorial clarity, is appreciated. — 
_ A word about the relationship of the elasti-plastic methods of this paper a 
‘ _ and a preceding paper 2) is inorder. The method of the preceding paper is 
~ analogous to the column analogy or elastic center methods used for elastic so- 
_ lutions, and it should be used for an elasti-plastic solution in any situation le “m ‘4 
_ which the column analogy is suitable for an elastic solution. The method of ie 
the present paper is analogous to the Cross moment distribution method used 
_ for, elastic solutions, and it is suitable for elasti- -plastic solutions in oil 
eS for which moment distribution is suitable for elastic solutions. = = © 
authors thank the discussers for their imaginative and valuable contri- 
butions to our of | on behavior, analysis, and design 


structures. 
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Closure Charles Bowman and Vaughn E. 


é. BOWMAN, 1 and VAUGHN E. HANSEN, 2 F. ASCE.—The 
- wish to express their thanks to Messrs. Sarpkaya and Macagno for their inter- a 
‘The purpose of the writers’ paper was to introduce a method of simplifying _ 
- dimensional analysis and to clarify some of the essential features. Before lJ 
dimensional analysis in any form can be efficiently applied the user must have _ q 

Dall knowledge and experience in the field of use. It can be a powerful — 


a 
pretations are made the effectiveness of the tool is s lost. 
z _ The vectorial analysis as presented by Sarpkaya is a ee | 
4 _ which allows for solution of problems which have four unknowns and three 
equations which can be readily seen and a fourth equation which is a 
‘ ¥ an 4 of two or more of the first three equations, and which cannot be readily writ- 
ten down. This type of problem occurs when (n-r) equals four, and there are 
= three non-repeating variables of which one is a combination of two of the four 
variables. This being demonstrated by the following fluids problem 
which involves the following variables: g, F, Q, L, V, and p. Qbeing 
= function of L and V. There are two dependent variables present and when using — 
eins proper procedure with repeating variables in the matrix ‘method with F, L, and . 
_s as the basic dimensions, it is difficult to solve for the four dimensionless a 
g 4 _— paxameters. . By using the vectorial analysis the solution becomes very easy. = 
: a B is combined with F, giving a resultant of LF, thus the polygon is complet- a 
by Q2. analysis gives the following numbers: Nr, NF, Ne, and 
when all of ihe variables are is a of Ne 


The value of (n-r Yr) is 
e merely a guide and will indicate the number of dimensionless parameters a 
+ which can normally be expected where many variables are involved. This - q 
ee equation may or may not be valid in all cases. « There is no way of identifying " 


the types of variables involved by this equation, therefore, it is an approxi- 


mation. If three variables are used, any one of the following, force, “Tength, f 
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mensioniess: parameter c can ‘be formed with three by of 
_ the variables and their exponents even though n-r equals zero. If only one of a 
- the variables contain a force, length, or time factor, and the rest do not, then _ 


n-r will equal zero nee a dimensionless eee cannot be formed from all > 
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Pd... of the method described in the paper are enhanced by the clarity of 
presentation and detailed discussion of its many possibilities. _ Of particular — 

_ interest is the rigorous treatment of convergence and stability conditions. 
— purpose of the present discussion is to point out a refinement of the 
_ method, which increases its accuracy without unduly complicating the ee j 

: cal work. The refinement is basedon the well known analogy between moments 
in beams and displacements of dynamic systems.\!) In a beam shears are ob- . 
oo tained by ee the loads, and bending moments by integrating the shear- ‘i 


When the has more than one component, the analogy can 
_ be set up with each component. And when the motion includes rotational com- 
ponents, and it is important | to consider rotational inertia, the acceleration, — 
velocity, and displacement in the analogy become angular quantities for those 
ec: the writer’s knowledge the most general and expedient procedure for ; 
aaa _ computing bending moments in a beam subjected to distributed load consists — 
_in dividing the beam into equal segments, approximating the load by means of 
a series of second degree parabolic arches, finding concentrated forces stati- 7 
cally equivalent to the parabolic loads, adding these to obtain shearing forces, 
and adding these in turn to find the bending moments. (2) _ Effecting the trans- 
bs. _ formation required by the analogy, the formulas used with equivalent concen- e 


1 Dir., Inst. of Eng., » Univ. of Mexico, Mexico, D, F, — a 
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n n+ 

¥, is the average velocity in the intervalt,,t, +h. 

Pe. expressions are used in the same manner as the e author’ Ss | Eqs. ; @) 


ae “acting at time to once the | nth displacements and velocities are known. If one i 
assumes computation of 7, (Eq. (2)), (Eq. (1)), and vn+1 (Eq. (3)) is 

The assumed acceleration is then revised. 

undamped the ‘may not be interest. ‘Eqs. (2) and (3), 


In may be some difficulty at n=0, for both aj and ag 
a have to be assumed. ‘This is overcome either by using the present. parabolic — ig a 

ie approximation asa refinement on a previously obtained solution by the author’s oe 
method, | or by using that method for the first interval and the present one for 
the rest, or by changing the computer’s program so as to assume ay 1 and _ 
— and revise their values untilconvergence is attained, proceding with assumption — a 
fa single acceleration in the rest of the time intervals in question. Te. ae cM 
A treatment of stability and rate of convergence similar to the author’ _ 
eads to his Eqs. (14) - (26) with Y = 1/2 and B = = 1/12. Accordingly the pres- ce 
“ent procedure has the ‘disadvantage > of permitting instability if the time inter- j a 
 Valis chosen only with a look to convergence. . The difficulty is overcome either 
using this method to refine solutions already obtained by the author’ or 
fea by making a separate check on stability, or, when possible, by resorting to rm 5 
revious experience with the type of problem athand. 
~The advantage of the present n method lies in its ‘Errors incurred 
in using the author’s procedure are of the order of h2. ‘They are reduced a "a 
é... order of h4 by the present variation. Hence the refinement is justified x, 2 
The pang of points required to describe accurately the ground motion = i 
so much so that even on modern computers the time it takes i : 


case, it may be preferable to use the more refined variant. 
i A word of caution is in order concerning cases of discontinuity in the ; ac- 


es in subscripts, or extrapolating the acceleration diagram as is done in simi- 
lar cases of discontinuity in beams. Otherwise the parabolic 
yields results less accurate than the unrefined method. 
. Frank Baron, “Successive corrections— pattern of thought, 
Methods of Analysis in Engineering, arranged and edited by Le E. Grintner_ 
} Macmillan New York, 1949), 127 
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M. Newmark, for computing deflections, ‘moments, 


_ and buckling loads,” Trans. ASCE, Vol. 108 (1943), pp. 1161-1234. rag par 
" T. P. Tung and N. M. Newmark, “Shears in a tall building subjected to. 
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ture? If not, is it satisfactory to base a design, not on i its expected perform- 


COMMENTARY ON PLASTIC DESIGN estan Dt STEEL: BASIC CONCEPTS# 


we However, as pointed out recently, the foundation of the plastic method of = 


design, being itself ‘subject to plastic” flow, appears to be on the verge of 


It ry thata satisfactorily = 
aden structure is expected throughout its useful life to carry “its intended © 7 
loads and otherwise performs its intended function”. This is no doubt a healthy _ x 
eed philosophy. However, is “collapse” an intended function of the designed struc- — 
_ ance, but on a condition — should almost never appear during its waekel 
“At working | load”, according to the Commentary, “the plastically designed — 
z structure is normally i in the so-called elastic range”. No matter what is meant 
~ by the so-calledelastic range, it is certain that such structures are “ normally” a 
hot in the so-called plastic range. . In other words, the plastically designed 
‘Structure is “normally” abnormal, and does | not behave e as the “plastician” in- 4 
tended it to. ‘This as asserts that the performance « ofa . structure is “normally” a 
_ divorced from the design in the plastic method. _ Can a design be judged to be 
if it takes very little account of its ‘normal performance? 
_ (2) There is a tendency in the Commentary ' to accept the “collapse ‘mechan- 
ism ” as the limit of usefulness. For instance, it is mentioned that ' Spore * 


n 
that the plastic method suits continuous beams and frames, which arethe major 


constituents of f any building, it is only fair that the occupier of such buildings Oy ; 
Ba should be consulted. | However, asa “plastician” admitted, 5 “the occupier of a 


building is of course more interested in its behavior under working loads than 


& 


| 
2 Su, H. L., “Design by quantitative of safety® "Tran, ASME, no, 59- 48, 


oe 3 Su, H. L., Disc. on plastic design, Peek. ICE, London, v. 13, p. 410, dul. 1959, 
H. L., Disc. on J, F. Baker’s 


§ sT 4, ‘Apr. 1959, pap. no. 2005. 
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tence of the Commentary, that is structure should be designed to perform 
its intended functions, and not to 
(3) The plasticity "developed in the actual steel structure mentioned 
i Section 1.2 of the Commentary obviously falls into a category completely dif- 
. ferent from the idealized, hence partly falsified, plastic behavior implied in a a 
_ collapse design. | This real plastic state will sooner or later disappear aft after S 
= yielding due to the effect of strain hardening as well as of relaxation. So, even 
Y in the Commentary, it is — that working load will only produce _— 
If this kind of real were to be incorporated in an 
-“plastician” would age rapidly and become “collapsian” eventually. Asa 
‘matter of fact, such plasticity has never been included in any plastic design by 


ian” 80 far. ‘It is hte that | steel must be subjected to some ‘kind = 


The plasticity mentioned i in the Secnienntuaned is in factan extensionof such pro 
cessing treatment. If such real plasticity can beused as an argument to justi 
can the plastic method, then the metal urgist can provide m many more reasons to 
_ (4) In comparison with many previous allegations made by plasticians, "the 
Commentary is quite modest in stating that for a large ‘majority of continuous: a 4 
and frames, a design based ‘upon a reasonable factor of safety (load “a 
factor) against this ultimate load provides a more appropriate owen than a 
ever, this Vague | statement to have been produced only to evade an es- 


One of the favorite arguments advanced | by * plasticians” ” is that theirs is 


load to the working load obtained by their method is “true” oy at least ' “truer” a a 
- that obtained in an elastic analysis. na 
_ It has been pointed out that whether a structure is safe anduseful cannot be : 
“judged on the basis of the collapse mechanism because a structure may have 
= needed serious repairs long before such a mechanism is developed. Further- ‘ 
ie more, corresponding to every pattern of loading, there will be a collapse | 
Ay mechanism. So long as the loading patternis | unknown, the mechanism — 
an uncertainty, therefore the load factor can hardly be “true” ; , nor more ap 
re. than the safety factor used in the conventional elastic design. = 
This also denies the necessity to detach the “allowable stress” concept from 
design, which « can be retitled the “allowable load” method without in- 
cn justice. It will be illustrated later that the difference between the elastic and _ 
i =e the plastic method of design is indeed no more than the difference between the se 
= word “stress” and “load”; and can never be polarised as “true” and “untrue”. _ 
6) The factor of safety has been | repeatedly and quite correctly called the fa 
“factor of ignorance in the past. No matter how cleverly disguised, the plastic 
method of design cannot escape from the reign of ignorance. Is there 
wy Treason why the load factor should be 1. 1.85 or whatever the value? - ‘The “plas- 
tician” can only reply, “past experience”! shows exactly the primitive 
_ The focus of this problem is directly and closely connected with the philoso- 
aber’ of une the principles of design. In many design problems, a great 


apse”. Moreover, this tendency of accepting “collapse” as the u 
2 
ae 
| 
ie 
aw | 
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- ‘Satisfy | the intended purpose. Generally speaking, this is a statistical prob- 
a ae lem. It happens that one of these problems (structural mechanics) can be pre- 
i see dicted with some certainty. Overwhelmed by this success, the engineer tends 
exaggerate the importance of the physical problem and at the same 


ignorance by introducing an artificial factor to safeguard the design as is well 
bes The problem of design is _s different from that of analysis . In both cases, 
- some predictions are made. In the problem of analysis, predictions are aN 
concerned with well-defined phenomena, and therefore can be made with accu- 
racy. In the problem of design predictions must be made on factors of various - a 
natures: physical, statistical, and economical. These fa factors may of them- 
‘selves be difficult to define, and their behavior may be so complicated that no- a 
_ rigid rule governing them can be found. Whether a structure is elastic or 
A _ plastic is only a, but not “the”, design problem. If the conventional elastic - on 


sign is unsatisfactory, it is far from fair to victimize elastic analysis 
{ ty discrimination. Many objections to the elastic method are in fact not physical _ 
statistical in nature. The has indeed expected too much and ex- 
acted: too ‘much from it. 
isa realistic method in of so ‘many abuses in ‘the past. 
(6) Many defects of the elastic method have been repeatedly used to ‘support 
e 


plastic method. This development seems rather illogical: the fact that 
White is a fool can be used justify that Mr. is a wise man, 
‘Mr. Green, Mr. Brown, etc. must also be wise men. a 
_ On the other hand, while attacking the weaknesses of me conventional elastic 
method, the “plastician” forgets that his collapse diagram cannot overcome 
such shortcomings too. The plastic method is merely a solution toa pie ll 
_ in structural mechanics, which only forms a facet of the design problem. It is 
oo that when a structure collapses, many effects and considerations become © * 
immaterial. However, these effects are real whereas “collapse” is only a con-— 
venient imagination. Is it possible to replace reality by sheer wishful thinking? 
_ Such real effects as the sinking of supports, non-rigid connections, etc. will | 4 
_ always exist no matter what design method is adopted, and canonly be reduced 
or eliminated by improving methods of fabrication and construction. On the 3 < 
_ other hand, following the advent of modern computers, most of these effects ‘= 
can now be incorporated in the elastic design; whereas such inclusion of real Hae 
effects in the method, possible, will make it very complicated and a 
eee (7) In the Commentary, it is : rather consoling to note that “plastic design ae 
a not a technique that is intended to replace all other design procedures”. — Un- “= q 
only fatigue and are mentioned as the others. seems 


and the plastic method be to designing structures 


ing, wire drawing, metal ete. to which the yield mechanism is indis- 


 pensible, the | theory indeed avery and 


a. 
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a . bution etc. are physical phenomena. If there exist parallel phenomena in s ea 

certain domain of th world, then th se of pia d on) absolutely 

* 
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The structural designer cannot. “enjoy the same privilege s since 


nothing his hand should flow or collapse. If such disasters are to 
avoided merely by means of an artificial factor of ignorance, the designer be- 
clouds his client as well as himself. It seems highly desirable that a much 


ought to to be ‘given now. That is, engineers should not beo over- 


— 


_  HSUAN-LOH SU. .1—(1) In Section 4 


ble for plastic design. this s sidesway is ; along» the “unconsidered third di- 

- mension, that is, the direction normal to the drawing, then whether plastic ie 
= is practicable depends upon how “multi” the number of redundancy 
_ andthe number of possible hinge positions are. For a structure of 24-redundancy © 
and 50 possible hinge positions, the total number of different combination of _ 
plastic mechanism will be 50!/25!25! + 1015, Is this problem still suitable for 
™ this sidesway is in the plane of the drawing paper, then these - walls and 
- bracings will not only prevent sidesway | but also handicap the formation of | 
_ plastic hinges. So the question of composite action emerges. Has the plastic — 


i method beendeveloped to this levelyet? Can it it ever be? If not, is | this problem a 
still suitable for plastic design? fas 
(2) The Commentary stated in Section 4.4 that at plastic hinge locations ~ 
punched holes etc. should not be permitted. This immediately raises the question oie 
“where”? The location of the plastic hinge is isa function of the 


a _ In the same Section, it was stated that “design details should be such that ze 
ae the material is as free to deform as possible”. This requirement seems a 
7 little strange, | because the opposite is more desirable. There is no pointin 
_ making an actual structure collapse so that the collapse design mayhave are- 
- alistic response. On the contrary, it is a designer’s duty to see that the de- : 
signed structure will be stiff enough and cannot collapse wire 
(8) Under the heading “loads”, section 4.7, the difference between elastic 
_ ne and plastic design is insufficiently explained. However, in both cases, a ej 5 
4 fi of spor gene is introduced. Consequently, both methods of design are unsatis- 
Saad (4) The Section dealing with load factors needs drastic revision because they 
i all factors of ignorance. A method for computing a realistic factor of a 
’ safety 7 ‘upon the theory of statistics and of probability has been +. acer 
recently.3 
2 Su, H. L, Discussion on the Commentary, paper 2091, Proc, mei aS 
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various factors the such as uncertainties in analysis, 
in workmanship, in stress concentration, etc. mentioned in the Commentary. Z 
sufficient statistical dataof these uncertainties and variations are available, 
computed safety factor will be very reliable with respect to a givenproba- 
bility of safety. If not, , the « designer c can still use his ; own judgment to — hes 


seems quite suitable for assessing “the value of a second hand car,” as one of 
the critics put it. The latter, in spite of making some use of statistics, bas 


the other by a a Committee of ASCE. The former is rather artistic 
been s started ona not very promising line. . Even the advocates of this: methot 


_ factor of safety will be ‘supplanted by the concept of probability of ot : 
(p. 1394, ASCE Trans. Vv. 121, 1956) and that “ a great deal of work remains ie 


> (p. 22. paper 1316, Proc. ASCE he 
 v83, n.ST4, Jul. 1957) It should be pointed out that such pessimistic views are 
= a consequence of the “margin of safety” concept. © The actual statist 
cal treatment of the safety factor, as has been shown3,4 is very simple and Ps 
ee Few criticisms can be made upon the details of Chapter 5, which is de 
voted to experimental verification of plastic theory. _ However, in concluding: 7 
this Chapter, it was stated in the Commentary that “in view of the notable — 
: agreement between plastic theory and the results of these tests, the applica- 
_ bility of the plastic method to structural design problems involving continuous 4 
steel beams and frames is demonstrated”! pos 
_ The notable agreement between plastic theory and the test ‘deeute is a fact 
. ey which can at most confirm the acceptability of the theoretical idealization of _ 
_ the plastic behavior of steel. In order that the applicability of plastic method fa 7 


_ to structural design can be “demonstrated”, there must also be some “notable i : 
= reality. ‘Unfortunately, there is no such “agreement” between the bare test © 
he frame, even the full-scale one, and the actual structural framework with brac- 
r ings, walls, and other infilling panels; nor between the ne loading pa — ened in _ 
the test and that likely to act on an actual building. 
‘There exists a link between a theory and an experime r as ‘the pure 
2 : physical problem is concerned, this is the only li link. The | «plastician” should 
be congratulated on his achievement in n defining and verifying one of such links. 
_ Nevertheless, there also exists a link between experiments and reality. This 
a second link is of a statisticalnature and can never be proved or demonstrated 3 
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SHRINKAGE, SWELLING AND CREEP IN en ae 


A. M. Neville and Keith Jones 


is. 


A. NEVILLE. .1—Little can be added to this expose but the 
. a might be interested, in connection with his observations on “Drying” to hear of 
some of the writer’s test results. _ Cylinders, 2 in. . in diameter and 9-1/4 in. > 
long, made of standard (British) mortar but with different Portland cements a 
were found to lose weight when stored continuously at a relative humidity of a 
95%. . They also exhibited shrinkage whose value between the ages of 1 and = 
‘months was | approximately 80 x 10-6, which corresponds to to a volumetric con- — 
4 ‘For the mix proportions used, the weight of water occupying this — 
: _- represent barely 0.01% of the total weight of the specimen, and yet tthe - 
_ actual loss of weight varied between 0.18% and 0.36%. It is known that the a ; 
hee magnitude of the actual shrinkage is affected by the restraining effect of the ‘ 
 non- shrinking aggregate and unhydrated cement grains, and by the action A a 


elastic forces at the Points of contact of solid particles, We can see that due 
_ to this loss of water or due to the self desiccation of the paste that exist in i 
Ss Be cement paste voids, which allow internal water movement without any apparent — 
: _ loss of water, and it is the transfer of water to and from the cement gel that is 
involved in reversible shrinkage andincreep. | 


ee may be of further interest to record that the loss of weight of loaded and — 
i. ou nloaded specimens was found to be sensibly the same. This was true for 
’ . stored ata relative humidity of 95% and also for specimens stored 
at 32%. In other words, any movement of water caused by creep can be in- 
ternal only. It may be interesting to recall that the author found the moduli of 4 
elasticity and Poisson’s ratios of moist and dry cements to be the same. ts 

4 ‘Thus, as mentioned before is likely to be related to internal 


JONES, 2 F. ASCE.—The author’ proposal of an “active water” theory 

for the explanation of shrinkage, swelling and creep in cement requires, as 
3 Mr. Hrennikoff recognizes, further research to determine parameters relating © ian 


_ water film thickness change to specimen length change. 
ol lt appears to this writer that the lack of agreement between the author’ : ee 
oa theory and his experimental results is due to failure to account for ; a large — 


percentage of the volume of “active” water. The need for a parameter 


oy 2 Engr., Office of Assist. Commissioner and Chf. eat, 1 os ‘Sie 1 of Reclamation, : 

Be Dept. of the Interior, Denver, Colo, 
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- treatment of the problem which is is based o on n spherical particles of sic size. 
heme of the 
 —_—_ of a cement specimen. The variation in sizes and ‘shapes of the 
= does not violate the concept of amorphousness and has been sub- 
Stantlated by electron microphotographs shown in the author’s reference (4). 
_ Assume that the specimen represented by the scheme of the ain 
ia figure of this discussion has beendried to the point of steepening of the volume- 
ie length change curve as described by the author’s Figs. 2 to 7. _ The author’s 
theory is that the remaining water is in films surrounding | each particle and 
this he has designated * “active” water. Suppose that drying of the specimen is > 
_ continued through the steep portion of the volume-length change curve driving — - 
off the active film water. Referring to the accompanying figure, only the re- me, 
_ duction of the films on the particles labeled “A” in contact and confined along ’ 
the length would contribute to the decrease in length of the specimen, while the v 
ay x diminishing of similar films on particles labeled “B” not closely confined and 
a free to move without affecting nearby particles would have no effect onshrink- 
_ age of the specimens. However, the films on the “B” particles would contribute i 
- the volume of water evaporated from the specimen, and the result would be 2 é 
ae a shrinkage smaller than that attributable to the reduced film thickness when 
- the author’s method of comparison is applied. The magnitude of the film 
_ decrease- -shrinkage ratio of that specimen as computed by the author’s method, ; 
_* depend upon the film volumes on the “B” particles as compared to those < 
‘Next, consider a specimen in whichall particles are incontact with at 
four other particles so that each is confined to one position with respect to me 
oa is equivalent to making all particles “A” particles. In that speci- 
men any active water which is evaporated would have an effect onthe shrinkage 
 € the specimen and w would correspond, except ‘for the variation in size and — 
shape of particles, to the structure assumed by Mr. Hrennikoff. _ <i 
ae condition necessary to produce such aspecimen is that a suitable gradation 


are p: present and so ) arranged that no void that could costhia the: smallest 

cle exists. A structural scheme approximating that condition or the condition — 
assumed by the author may have predominated i in the specimens made of cement i a a 
_ finer than #325 mesh which failed to survive the tests. Specimens ‘made of such - = 
"graduated cement particles should be less permeable to moisture and agreater Hs 
differentialin swelling during wetting and shrinkage during drying would neces-_ 
«sarily exist between the core a and the surface. If too rapid wetting or drying 
were permitted, early failure of the could’ be expected. ‘The film 

_ decrease- ~shrinkage ratio, then, should be related to the. gradation in size of | 


It is, therefore, this writer’s opinion tt that the active water thatdoes notcon- 


tribute to overall ‘swelling and shrinkage « of the specimen can be accounted for 

in films on particles which are freeor partially free to move in poresor cavi- 
‘ties and which are of irregular shapes so that more water is in a film ona aie 
particle than would be in a film ona spherical particle of the same volume. 
_ Although this writer cannot offer a quantitative parameter expressing the — 

5 effect of active water as related to the gradation of particle sizes, the experi- | 


- ments and theory offered | by Prof. Hrennikoff would seem to demand that some 
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Soyo DIVISION ACTIVITIES 
ENGINEERING MECHANICS DIVISION 
= of the American Society of Civil Engineers 


7 = As the years _roll by the Engineering Mechanics Division continues . 

: healthy growth. The number of papers submitted to the Division for publica- 


—1955- 56 


=) 
This increase in the number of the to in- 
augurate a new service to Society members; i.e., to the scheduling of a special ; 
spring meeting, described elsewhere in this newsletter, at which about a dozen = 
of these papers will be presented. Since it is to be centrally located at Purdue 
and held in conjunction with a student ASCE Conference, it is to be hoped that — os 
of you in academic life will consider bringing some your students 


= The “Commentary on Plastic Design in Steel” prepared by the Joint WRC- 
Committee on Plasticity Related to Design is being delivered on sched- 
a ule. Progress Reports 1-4 were published in July and October. Parts 5-7 are i 
scheduled for the e January and April issues of the EMD Journal res 
The Fluid Dynamics Committee preparing study on “Stratified Flow” 
“ and a new task force is being organized to consider “Flow Past Objects Im- rae 
. 3 Similarly, the Committee on Structural Dynamics is working on a manual < 


this field and should have it available soon. 
ae all some 75 members of the Division are working on its ait and 


for and criticisms will be welcome. 


Chairman, Executive Committee a 


_ Note: No. 1960-2 is part of the copyrighted Journal of the Engineering Mechanics Di- 
Proceedings of the American Society of Civil Engineers, Vol. 86, EM 


Coypright 19 4009 by the American Society of Engineers. 
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Te Engineering Mechanics Division of ASCE is planning to hold a Confer- 

ence on Structural Mechanics at Purdue University, Lafayette, Indiana, on 
% _ Thursday and Friday, May 5 and 6, 1960. This is the first conference to be 
4 sponsored by the Engineering Mechanics Division. Papers are being selected a 
4 from those recently submitted and already approved or currently being re- 
viewed by appropriate committees of the Division. Professor Bruce 
University of Michigan, Program Chairman ' for 1 1960 for the Engineering Me- 
Me 7 chanics Division, is acting for the Division in the. scheduling of the papers. mo 
The technical program will consist of three half-day sessions, one during 
= afternoon of May 5, and two on May 6. The morning of May 5 will be re 
ey aside for registration and inspection trips to laboratories on the campus. = 4 
a joint banquet with the Student Chapters of Notre Dame University, Rose a? 
- Polytechnic Institute, and Purdue University, and the Indiana Section of ASCE > eS 
will be held at the Purdue Memorial Union on the evening of May 5. oe Paz 


_ Programs will be published and distributed to members of the the ieee 


New in the field of strechesal s of ‘computers have 


_ continued to take place since the first Conference in Electronic Ceppeeten, 
_held in Kansas City, Missouri, in November, ee diy, 
_ Planning for the Second Conference is underway. — This conference is 
a nana for Thursday and Friday 8-9 September, 1960, at the new Pittsburgh | =. 
Hilton Hotel. Conference sponsors will be the Structural Division 
on Electronic Computation and the Pittsburgh Section, A.S. 


Four main themes are planned: 


including such areas as optimization of of s structures and least- ast-cost design. 


> 


5 ee” Design Applications - Outstanding applications in the field of design, — 
* 


2. Structural Analysis Emphasis will be | placed on n methods 
adapted specifically for use oncomputers. 
Mathematical Methods - New and techniques in numeri- 
‘analysis | methods for ‘el 
"Professional Problems - Education tion of ineers in the 


inquiries regarding the should be addressed to 
Chairman, Program Committee. Nathan M. Newmari Illinois 


_ WEST EST COAST CONFERENCE 
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jj. ‘ina preliminary announcement, Professor Gere, Secretary, West Coast § 
Committee, announced that the 1960 West Coast Conference of Applied Me- _ 


will held at California of Technology on Jun 

: This is an annual meeting sponsored jointly by the West Coast Committees, 
Engineering Mechanics Divisions of and ASME. 


=o 


Second on Structural Mechanics will be held April 5 
a 6, and 7, 1960 at Brown University, Providence, R. L, under the joint spon- 


_Sorship of the Office of Naval Research, Department of the Navy, and of Brown P ; 


t University. The Symposium will be devoted exclusively to the field of plas- — 
a ticity. The program will consist of critical surveys in selected areas and of 
_ reports on ‘original research, with ample time for discussion. The organizing © 
a ‘consists of Professors E. H. Lee and P. S. Symonds (co- chairmen), — : 
_ Professors D. Drucker and W. Prager. 


EW PUBLICATION POLICIES 


organization a number of changes regarding publication policies and proce- 
dures. The proposal appears to have been influential in the changes ¥ 
at the past annual convention. _ The sense of the changes is summarized in the 
following letter recently mailed Division Chairman by William Wisely, Ex- 
At the Annual Convention of Society in Washington, D.C., October 19-20, 


an 2. Accepting only the papers that can be recommended for inclusion in the 


3. Assembling, in the annual bound volume of “ Proceedings Papers” 
“Transactions”), all Journal papers that have attracted constructive discus- 
a iad implement this new program will require the active cooperation of ae ; 
_ members, Division officers, and authors, and the purpose of this letter is “a ji 
= your objective support. — It is anticipated that there may be some — 


“for ‘clarification in the details of the pesca and these will be issued 


"Publication 


“Discussion Rules* 


use of your Publications Committee Chairman and his reviewers. “Subject 


minor adjustments, the new procedures give advancing 


— 
prove the content and appearance of the Division Journals by--- 
— 


still: further, the high | status that “the ‘Society’ s have enjoyed 
more than 100 years. They can — with your continued encouragement 


AL REPOR a 7 
following report of activities for the Engineering Division 
- October 1, 1958 - September 30, 1959, with appended committee reports, was es 

_ submitted at the Washington, D.C. meeting by the retiring chairman: | Ridin 

“a This year the Engineering Mechanics Division welcomed a _ vigorous new 
_ Contact Member, C. W. Britzius, and Secretary, E. Wenk, Jr. and gave its 

heartfelt thanks to retiring Contact Member, R. Sherlock and Secretary, Merit | 


‘Then the Journal of the Mechanics Division and the Journals of 
. the other Divisions with a strict and careful review procedure will be able to 
be their appropriate place on the national and international scene. a ; 
a Executive Committee met in New York on October 14, 1958 and again . 
a on April 3, 1959 in Ann Arbor. _ Gradually the Division becomes more like the © 
aa in its administrative set-up. We now have an active Newsletter under 
: the able editorship of Professor Donald L. Dean of the University of Kansas. 
ce _A Research Committee is functioning under the chairmanship of Dean John S. a 
i“ McNown. Publication and Program Committees are well established and plans 


_ for the future are far advanced. One of the Technical Committees has a tack . 


The reports of the Committees which are give clear of 
‘most of the work of the Division. Among items not included is the matter of 


separation into Technical Committees. ‘Fluid mechanics in particular seems 


slighted both organizationally and b; by prospective au authors. Meteorology might 


The Division co-sponsored the West t Coast Conference of Applied Mechanics 
: a the Applied Mechanics Division, ASME, and the Midwestern Conference | 


_ on Fluid and Solid Mechanics with a number of interested groups. = “- 
A Technical Procedure Conference was not held this year, but I attended — 


and attitudes showed considerable progress along lines we have advocated for ays 2 
80 many years from our noisy | but narrow Sen pe as the Committee on | Ap- a 
_ Each retiring Chairman has pointed out that the major - share of { the work is 


Division. 


~ 
i 
ivities and publica- 
importance of the Journals over all the other va 
Mu _tions, Headquarters staff is doing an excellent job un ee 
— 
— 
— 
— @ 
| 
r r effectiveness and desir- 
April 17, 1959. Much discussion centered 
abili af the Contes Memher stom Publications 2 
the Committees and most of all by the Secretary. They, it 
pss their gratitude to the many who helped them — 


Engineering } Mechanics Division 2--5 


Merit P. White (1963) Incoming Member ae 

Wenk, Jr., Secretary 
C. W. Britzius, Contact 
Daniel C. Drucker (19589) | Chairman 


s. Archer 


Fluid 


Mechanical Progerties ki G. Dietz = 


os tee with Dan H. Pletta as Chairman and the Program Committee with Egor P 


Popov as Chairman for 1959 and Bruce G. Johnston as Chairman for 1960. 


J 


Board 


P. White 


"Submitted 
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